Existence and stability of stationary solutions to the full 
compressible Navier-Stokes-Korteweg system 
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Abstract 

This paper is concerned with the existence, uniqueness and nonhnear stability of stationary 
solutions to the Cauchy problem of the full compressible Navier-Stokes-Korteweg system effected 
by external force of general form in M.^. Based on the weighted-L^ method and some elaborate L°° 
estimates of solutions to the linearized problem, the existence and uniqueness of stationary solution 
are obtained by the contraction mapping principle. The proof of the stability result is given by 
an elementary energy method and relies on some intrinsic properties of the full compressible 
Navier-Stokes-Korteweg system. 
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1 Introduction 

In this paper, we are interested in the following nonisothermal compressible fluid models of Korteweg 
type, which can be derived from a Cahn-Hilliard like free energy( see the pioneering work by Dunn 
and Serrin [2], and also [3 HE]). 

{pt + V-{pv) = G{x), 

{pv)t + V-{pv^v) = V-{S + K)+ pF{x), 



P[e + - 



+ V- 



pv 



v^ 

^ + T 



(1.1) 



V • {aVe) + V ■{{S + K)-v) + pv F{x) + H{x), 



Here (x,t) € M^ x M+, p > 0,^ = (^1,^2,^3))^ > and e denote the density, the velocity, the 
internal energy and the temperature of the fluids respectively, a is the heat conduction coefficient. 
F{x) = (Fi(x), F2{x) , F^{x)) , G{x), H{x) are the given external force, mass source and energy source, 
respectively. The viscous stress tensor S and the Korteweg stress tensor K are given by 



Sij = {p'V ■ V - P{p, e))6ij + 2pd 



%]\ 



K, 



«j 



'-{Ap"^ - \Vp\'^)6ij - Kdipdjp, 



(1.2) 
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where dij{v) = {diVj + djVi)/2 is the strain tensor, P is the pressure, /i and /i' are the viscosity 
coefficients, and k is the capillary coefficient. Notice that when At = 0, system (jl.ip is reduced to the 
compressible Navier-Stokes system. In this paper, we consider the case of e = CyO, where Cy is the 
heat capacity at the constant volume. Our basic assumptions are as follows: p, 6, k, /x, p' and a are 
the constants satisfying p, 6,K,n,a > and |^ + /i' > 0; Cy > is a constant and P = P{p, 0) > 
is a smooth function of p, 6* > satisfying Pp{p, 6), Pe{p, 0) > 0. 

The main purpose of this manuscript is to study the nonlinear stability of stationary solutions 
to the Cauchy problem of the compressible Navier-Stokes-Korteweg system (1.1). It is convenient to 
study the Cauchy problem for the following form which is equivalent to (jl.ip for classical solutions, 

(pt + V-{pv) = G{x), 
pvt + p{v ■ V)v = fiAv + (m + /i')V(V • v) - VP{p, 9) + KpVAp + pF{x) - vG{x), /^ gx 

pCy {Ot + {v ■ V)e) + ePe{p, e)V ■v = aM + ^{v) + $(/9, v) + H{x) + y G(x) - CsjG{x)e, 

with the initial date 

{p,v,e){t,x)\t=Q = {po,vo,eQ){x) ^ {p,Q,6) as|x|^+oo. (1.4) 

Here 



■^{v) = p'{V ■ v)^ + 2pOv : Bv, Bv = idijiv))lj^^, 
$(/), v) = K (^^ + pAp) V • v - K(V/9(g) Vp) : Vv 



(1.5) 



The stationary problem corresponding to the initial value problem (|1.3|) . (|1.4|) is 
fV-(H = G(x), 

(v ■V)v = - \pAv + ([1 + ^')V(V • v) - VP(p, e)] + kVAp + F(x) - -G(x), 
P P 

(v ■ V)e + ^^MMl V . ^ = J- jaA^ + ^{v) + Mp, v) + Hix) + ^G(x) - Gs^G{x)e } , 
pUv pGy [ 2 



(1.6) 



Before stating our main results, we explain some notations as follows, which are borrowed from 
[6] and [7]. 

Notations: Throughout this paper, we use the standard notation in vector analysis. For example, 
we put for scalar u, vectors v = {vi,V2,V3),w = {wi,W2,W3) and matrix / = {fij)i<i,j<3- 

A v^ ^^^ A /A A A \ / ^\ v^ du 

Au = 2^— ^, Aw = (Awi, At;2, Avs), [v \/)u = 2^Vi-—, 

i=l ^^i i=i "^^ 

{v ■ V)w = {{v ■ S/)wi, {v ■ V)w2, {v ■ V)w^) , 
V^n = {9>||a[ = k], V^v = {d'^v^\\a\ = A:,i = 1, 2,3} , 

^ dxi ' I ■^ dxi ' ■^ dxi ' ■^ 9x,- 

i=i \i=i ^ j=i ^ i=i J 
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Here a = {ai,a2,as) is a multi-index, |a| =01 + 02+03 and d^ = 8'°^' / dx"^ 8x2^ dx'^'' . Moreover, 
we also use the notations 



^2 



^ , du du du\ f d 

vu = - — , - — , - — , Vn 



u 



dxi dx2 dxsj \dxjdxiy ^^^.^^ 

and denote ^ by diu or Ux^ without any confusion. 

Next, we introduce some function spaces. Let L^ denote the usual L^ space, put for scalars ui,U2 
and vectors v = (t'i,'V2,- • ■ ,Vn), w = {wi,W2, ■ ■ ■ ,Wn), 

1 - .1 



\ui{x)\Pdx] , II^'IIlp = ^ ||fj|lLp ' (l<p<oo) 



\Ul\\LP 

||'Ui||l°o =sup|ui(x)I, ||w||l°° = max ||fj(x)||Loo, (^1,1*2)= / uiU2dx, 

K3 l<i<n 7r3 

/ \5 

{v,w) = '^{vi,Wi), \\v\\k={ ^ I|V'-u|m ^i'^*^ II • II = II • IIl2 

1=1 \0<l<k J 

H'' = {ue Ljj \\u\\k < 00}, H'' = {u£ L}J Vn G H'''^}, 
where u is either a vector or scalar. Further we put 

W'^^l = {{a,v,^)\a€ W, v G i^^19 G H^}, 

and 

||(o-,w)|| = ||cj||fc + ||v||«, \\{a,v,'d)\\j^k,l = \\cr\\j + ||f lU + \W\\l- 

Definition 1.1. 

I^ = {ae H''\ \\a\\jk < e}, J^ = {v e H^\ \\v\\jk < e}, N^ = {?? G H''\ W^j^t < e}, 
where 

k 

Ikll/^ = lkllL6 +E 11(1 + |a^ir(VV, V^+V, V^+MII + 11(1 + |x|)V||^^ + 11(1 + \x\fVa\\^^ , 
i/=i 

k k 

\\v\\j^ = Mp + E 11(1 + i^ir 'v'^HI . ii^iu^ = ii^iii^ + E 11(1 + i^ir''(v'^^' '^''^'^)\\ - 
i/=i i/=i 

and II • lljfe is defined by 

1 

hWp = \H\l^ + E 11(1 + I^D'^'^v'^ILoc + 11(1 + kl)'v 



'"IL- 



u=0 
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Moreover, we put 

Af '' = {(a,t;,^)|a G Ii,v € Jl^ e 7Vi,||(a,7;,T?)|U„M < e}, 

Ap-' = {(fj, TJ, T?) e Ap''| V • v = V • Fi + ^2 /or some ^1,1/3 

suchthat ||(1 + |x|)3yi||^^ + ||(1 + |a;[)-iy2||^i < e} , 
£ = {C/If/ = S/ ■ Ui + U2 for some Ui, U2 andsatisfies \\U\\c < 00}, 



where 



Mc = 5Z 11(1 + I^D'^'^V'^C/II + 11(1 + |x|)3([7, V[/)||^^ + 11(1 + \x\fUi\\^^ + ||C/2||ii. 



1^=1 



In this paper, we consider the case where the mass source G, the external force F and energy 
source H are given by the following form 





where Fi = iFi^ijix))i<ij<3, -F2 = (i^2,i(a:))i<i<3; Gi = (Gi,i(x))i<j<3, G2 = ^2(2;); Hi = 
{Hi^i{x))i<i<3, H2 = H2{x). 

Now we begin to state our main results. As [7], regarding p as a smooth function {P,6), Our 
first Theorem is concerning the existence of stationary solution to ()1.6p , and its weighted-L^ and L°° 
estimates. 

Theorem 1.1. Let p, 9 he any positive constants, and set P = P{p,9). There exists small constants 
Co > and eo > depending on p and 9, such that if {G, F, H) G T-L'^'^'^ and satisfies the estimate: 

\\{G,F,H)\\c + 11(1 + \x\)^V\G,H)\\ + 11(1 + |x|)-^G||^, < coe 

for some positive constant e < eo, then /il.6\) admits a solution of the form: (P, v, 9) = {P + a, v,9 + '&) 
where {a,v,9) E A^' ' . Furthermore the solution is unique in the following sense: if there is another 
solution {P + ai,vi,9 + i}i) satisfying 111.6]) with the same {G,F,H), and \\{ai,vi,'di)\\j^4,5,5 < e, then 
{ai,vi,'&i) = {a,v,'&). 



Next, we consider the stability of the stationary solution of (|1.6p with respect to the initial 
disturbance. Let {p* ,v* ,1!}*) be the stationary solution obtained in Theorem 1.1, then the stability 
of (/>*, w*,"!?*) means the solvability of the non-stationary problem ()1.3p . ()1.4p . Let us introduce first 
the class of functions which solutions of (|1.3p . (|1.4p belong to. 



Definition 1.2. 



0(0, T;W 



i,k,l\ 



{a,w,'d){t,x) 



ait,x)eG^iO,T;W)f]C\0,T;W-^), 

w{t,x) E G^{0,T;H'')f]C^{0,T;H''-^), 

^{t, x) e C0(0, T; H^) n C\0, T; H^-^) 
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Then, we have the following Theorem. 

Theorem 1.2. There exist C > and 5 > such that if ||(po — P* ^vq — v* ,6q — 'i?*)||4,3,3 < 5, then 
the Cauchy problem lil.3\) . |j.^[ j admits a unique solution (p, v, 9) = {p* + a,v* + w, 9* + "i?) globally 
in time, where {a,w,'d) G C{0,oo;'H'^'^^^), Va G L'^{0,oo;H'^), Vw,V^ G L'^{0,oo;H^). Moreover, 
the solution {a, w, f?) satisfies the estimate: 

\\{a,w,m)\\h,3 + f l|V(a,tx;,^)(5)||2 ds < C\\{a,w,^)ml,^,. (1.7) 

Jo 

for any t > and 

\\{a,v,'d){t)\\L°° ^ ast ^oo. 

The compressible Navier-Stokes-Kortewg system has been attracted many attentions due to its 
applications in fluid mechanics as well as mathematical challenge. A lot of mathematical results on 
such system have been obtained. More precisely, Hattori and Li |121 [T3] proved the local existence 
and the global existence of smooth solutions for the compressible fluid models of Korteweg type 
in Sobolev space. Danchin and Desjardins [TT] proved existence and uniqueness results of suitably 
smooth solutions for the compressible fluid models of Korteweg type in critical Besov space. Bresch, 
Desjardins and Lin [8j showed the global existence of weak solution to the compressible fluid models 
of Korteweg type, then Haspot improved their results in [U]. The local existence of strong solutions 
for the compressible fluid models of Korteweg type was proved by M. Kotschote [13]. Wang and 
Tan [15] established the optimal L^ decay rates of global smooth solutions for the compressible fluid 
models of Korteweg type without external force. Recently, Li [T7] discussed the global existence 
of smooth solution to the following Cauchy problem of the isothermal compressible fluid models of 
Korteweg type with potential external force. 

Pt + V ■ {pu) = 0, 

ipu)t + V-{pu^u) = V-{S + K) + pF{x), (1.8) 

_ {p,u)\t=o = {po,uo). 

Here F(x) = —'V<J3 with <j) being a scalar function and S,K are defined as in (jl.2p . He proved that 
there exists a unique stationary solution (p(x),0) to problem (jl.Sp if 4>{x) satisfies some smallness 
condition in the H^ norm. The nonlinear stability of the stationary solution (/5(x), 0) and the optimal 
L^-decay rate of smooth solutions to (jl.Sp were also proved in p!7]. Motivated by the work Y. Shibata 
and K. Tanaka |6] for the study of compressible Navier-Stokes equations, when the external force is 
given by the general form F = V • Fi + F2 and also mass source G appears, it is expect that the 
stationary solution is nontrivial in general. On the other hand, all the above results are concerning 
about the isothermal compressible fluid models of Korteweg type, for the nonisothermal compressible 
fluid models of Korteweg type, fewer results have been obtained. To our knowledge, the only available 
result for the nonisothermal case is [TU] , where the existence and uniqueness of strong solutions was 
proved in critical space. Based on these observations, we consider in this paper the nonlinear stability 
of stationary solutions to the full compressible Navier-Stokes-Korteweg system (jl.ip . 

Now we outline the main ideas used in proving our main results. The proof of Theorem 1.1 is 
motivated by the method developed by Y. Shibata and K. Tanaka [6j. Firstly, as mentioned before, 
we choose (P, v, 9) as the independent variables and regarding p as a smooth function of {P, 9). Then 
in the same sprit as [U], we need to establish the corresponding linear theory in the L^-framework 
for ()1.6p by employing the Banach closed range theorem. Compared with the case of compressible 
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Navier-Stokes system, the appearance of the third order terms VAo" and VA-!? in the velocity equation 
(2.7)2 result in more difficulties when we estimate the Li^ norm of the solutions to the approximate 
problem. In particular, an additional term V(V • v) appears in the energy estimate. To close the 
LP' energy type estimate, we frequently use the structures of the approximate system. Then by 
choosing some suitably space-weights and multipliers, the weighted-L^ estimate of solutions to the 
linearized problem is also obtained. In order to deal with the nonlinear problem, we have to derive 
the weighted-L°° estimates for solutions (cr, f,"!?) to the linearized equation ()2.8ip . The weighted- 
L°° estimates for v and •& can be deduced in the same way as that of compressible Navier-Stokes 
equations. However, for the weighted-L°° estimates of cr, we need to perform some delicate estimates 
related to the Bessel potential(see (j2.98p for detail). Moreover, the highly nonlinear terms ^{v) 
and ^{p,v) in ()2.82p are overcome by some delicate analysis. Having obtained the weighted-L^ and 
weighted-L°° estimates of solutions to the linearized problem. Theorem 1.1 follows by the contraction 
mapping principle. As for the nonlinear stability of the stationary solution obtained above, the key 
step is to deduce some certain a priori estimates for solutions to the initial value problem ()3.ip . ()3.2p 
in the H^ framework. Based on the properties we obtained on the stationary solution and some 
delicate estimates, we can deduce the desired a priori estimates. It is worth to point out that, for 
the compressible Navier-Stokes-Korteweg system (jl.2p . the appearance of the Korteweg tensor pVAp 
results in more regularity for the density than the velocity and internal energy (see (|1.7p ). In fact, we 
frequently use integration by parts and the equation (j3.ip 1 when we estimate the the terms containing 
VAo". As a result, the Korteweg term is split into the energy and the terms small in the L^ norm. 

Another interesting problem is to investigate the convergence rate of the non-stationary solutions 
to the stationary solutions when the time goes to infinity. As mentioned before, this problem has 
been studied by some authors for the isothermal compressible Navier-Stokes-Korteweg system with 
G = 0,F = — V(/) or without any external force(cf.jl7j. |15tll6j). But to obtain the convergence rate 
in our case, it appear to be more delicate since the stationary solution is nontrivial generally. We 
will consider this problem in a forthcoming paper. 

Before concluding this section, we also mention that the nonlinear stability of stationary solution 
for the compressible Navier-Stokes system has been studied by many authors. For the non-isentropic 
case, we refer to [201 EI] for the stability of constant state {p,0,9) in M^, [22j for the stability of 
nontrivial stationary solution {p*{x),0,6) in an exterior domain of M'^ and [TJ [23] for the stability 
of generally nontrivial stationary solution {p* (x) , v* (x) , 6* (x)) in M.^ and an exterior domain of M^, 
respectively. For the isentropic case, the interesting readers are referred to [HI [^ ESj for the stability 
of generally nontrivial stationary solution (p* (x) , v* (x)) in M^ or an exterior domain of M^ and [23] 
for the stability of nontrivial stationary solution {p*{x), 0) in an exterior domain of M^. 

The rest of this paper is organized as follows. In Section 2, we study the stationary problem. The 
non-stationary problem will be studied in Section 3. 



2 Stationary problem 

This section is devoted to the stationary problem (|1.6p . Take any two constants p,9 > 0. As 
mentioned in Section 1, by regarding p as the function of {P,0), changing the variables {P,v,9) -^ 
{P+a, V, O+i}), and rewriting the third equation by using the first one, ()1.6p can be then reformulated 
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as 



V ■v + —(v-V)a 
P 



P P 



where 



-fiAv - {fj, + /i')V(V -v) + V(T - K7iVAcr - k72VAt? = -p{v ■ V)v + /, 

-aA^ = -r]i{v ■ V)t? + mi^ ■ V)a + ^(u) + $ - r/gG + i^ + — G(x) - Cv(?? + 9)G, 



71 = PPP, 72 = We, PP = Pp{P,0), pe = peiP,0), 



(2.1) 



Vi = pCs, 



V3 



Ope 



(2.2) 



, V2 = — ._ 

PPP P PPP 

f = up (V(T • V'^pp + Vpp • V^o- + WppAa) + k {ppp - ppp) VAa 
+KP (Vi? • VVe + Vpe • V¥ + VpeAi?) + k {ppe - ppe) VAt?, 

$ = K [\\ppVa + peVt?p + p(Vpp • Va + ppAcr + V pe ■ Vi? + peAi?)] V • v 
-K [{ppVa + pg\7'&) ®{ppVa + peV^)] : Vv. 



Our goal of this section is to prove Theorem 1.1 by application of weighted L^-method to the 
linearized problem for (j2.ip . 



2.1 Weighted Li^ theory for hnearized problem 

We shall consider the linearized equation of ()2.ip : 

' V -v + ia- V)(T = g, 
-pAv - (/i + p')V{V ■v) + Va - KjiS/Aa - K72VA1? = /, 
-dA'd = h, 



(2.3) 



where a = {ai{x),a2{x),a^{x)), {g,f,h) G Ti '■^''^ are given. Throughout this subsection, we put 

/ = -(6i-V)ci + /, h = -{b2-V)c2 + h. 



and assume that 






(2.4) 



gj,h] €?^^'^'^ 6i,62,ciG J^C2€7V^ 



(2.5) 



:i + \x\){g, h)\\ + ^(1 + IxirV^ig, h)\\ + ^(1 + [xir+^V^if, h)\\<<^ 



(2.6) 



u=l 



u=0 



Zhengzheng Chen and Huijiang Zhao 



2.1.1 Solution to approximate problem 

First, we solve the approximate problem: 

V ■ V + (a ■ V)(T — eAo" + ea = g, 

-fj.Av - {fj. + Ai')V(V • u) + V<T - «;7iVAcr - k72VAi? + eu = /, (2.7) 

-aA^ + ed = h, 

in Ti"^'"^''^. In the following lemma, we prove some fundamental a priori estimate needed later. 



Lemma 2.1. Suppose that {a,v,'d) € H ' ' is a solution to \2.1^ . Then there exists two posit 



ive 



constants 6q = 6Q{'yi,'j2-,K, fi, fi',d) and eo = eo{'yi,'j2,K, fi, fj,',d) < 1 such that if 6 in ^2.4\) satisfies 
5 < 5q and < e < eo, we have the following estimate: 

|lV(a,z;,^)|||i^2 + e|l(^,^>^)ll2,i,i<Ce-^ll(5,/,/i)f + C|lV(5,/i)f. (2.8) 

Here, C > is a constant depending only on 71,72, k, ^, fi' and d. 

Proof. The proof consists of four steps. 

Step 1. Taking the L^ inner product with a and v on (j2.7P i . (|2.7p 9 , respectively, using integration 
by parts and canceling the term (Vo", v) by adding the two resultant equations together, we have 



/i||Vv||2 + {fi + ^')\\V ■v\\^ + e\\{a,v)\\lo 

= {9, <y) + (/, v) + K7i(VAa, v) + K72(VAi?, v) - {{a ■ V)(t, a). 

Differentiating (|2.7P i and (J2.7p 9. and employing the same argument, we have 

^||V2„f + (^ + /.')||V(V . vW + 6||V(a, v)\\% 



(2.9) 



(2.10) 



= (V5, Va) + (V/, Vv) + K7i(V(VAa), Vv) + k72(V(VAi9), Vf ) - (V((a • V)ct), V(t). 

Adding (fZTOD to ([IJ]) yields 

^Ji\\Vv\\l + {^l + ^l')\\{V■v)\\l + e\\{a,v)\\l^ 
1 
= Y. {^^''^' ^''^^ + ^^''^' ^""^^ + K7i(V'^(VAa), V^t-) (2-11) 

+k72(V^(VAt?), V^f;) - (^^((a • V)a), V^a)} = /^ + /2 + Ig + 14 + h. 
It follows from the Cauchy inequality that 

h < |lkf + a\\gf + 7]\\Aaf + CJgf, 

h <\\\vf + ail/f + 7?||V(V . ^)f + C7,||/f , (2.12) 

h < ??(||Vaf + llVAaf ) + C^||V(V • v)f, 

and 

h < r?(||V^f + llVA^f) + C^||V(V • v)\f. (2.13) 
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Here and hereafter, rj > denotes a sufficiently small constant and C^,Cr^ denote some positive 
constants depending only on e and r/, respectively. Moreover, the Cauchy-Schwartz inequality and 
the Hardy inequality imply that 



/5<C'(|((a-V)a,a)| + |((a-V)a,Aa)|) 



<C||(l + |2;|)a||Loo ||V<t|| i-r + t^ ||Act| 

< C(5||Va||f . 

Combining ([2lT]l - (f2Ji]) . we obtain 

^i\\Vv\\l + {^Ji + ^l')\\y■v\\l + e\\{a,v)\\l^ 

< Cr?||V(a,^)||i + C6\\Va\\l + C,||V(V • v)^ + {C, + C.Md, fW- 

Step 2. Differentiating (|'2.7P i . we get 

V(V • v) = -V((a • V)(t) + eVAcT - eVcr + Vg. 
which together with the sobolev inequality imply that 

||V(V • v)f < C{\\{a, Va)||Loo \\Va\\l + e^\\VAaf + e^\\Vaf + \\Vgf] 

< C {6^Va\\l + e2||VAa||2 + e^Vaf + \\Vg\\^) 



(2.14) 



(2.15) 



(2.16) 



Step 3. Taking the L inner product with Vo" on (J2.7P ';. we have from the Cauchy inequality 



that 

||Vo-|p + K7i||Ao-|p 

= fi{Av, Va) + {n + ^')(V(V • v),\/a) + K-f2{VA^, Va) - e{v, Va) + (/, Va) 

< -||Vct||2 + C {\\Avf + ||V(V • v)f + ||VAt?||2 + e'^Wvf + \\ff) 

Consequently, 

\\Vaf + IIAfjf < C {\\Avf + llVA^f + e^\\vf + \\ff) . 

On the other hand, it follows from (|2.7|) 9 that 

\\VAaf <C{\\Avf + \\VA^^ + \\Vaf + e^\\vf + \\ff) . 
Therefore, we have from a linear combination of (|2.18|) and (|2.19|) that 



|Va||^ < C {\\Avf + \\S/Ai}f + e^Wvf + Wff) . 



(2.17) 

(2.18) 
(2.19) 

(2.20) 



Step 4. By using the same argument as ()2.9p and ()2.10p . one can get from ()2.7p -? that 



a||Vi?||? + e 



(/i, t?) + (V/i, Vt?) 



< - 

- 2 



'—' + ^\\v'n' + ic,+c)\\hf 



which implies 



a||V^||? + 6||^||?<(a + C)|| 



(2.21) 



(2.22) 
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On the other hand, 

a||VAt?||2 = II - eV^ + Vhf < e^HV^f + \\Vhf 

Combining H^TU^i and (I223D, we obtain 



(2.23) 



(2.24) 



Thus, by some suitably hnear combinations of (|2.15p . (j2.16p . (j2.20p and (|2.24p and using the smallness 
of €,7] and 6, we can get ()2.8p . This completes the proof of Lemma 2.1. 

Now, we employ the closed range theorem to prove the existence of solution to (j2.7p . We introduce 
the operator A defined on D{A) C L^ into H^ x L'^ x H^ by 

A{a,v,^) = {Ai{a,v,^),A2{a,v,^),As{a,v,^)) 

where D{A) = T^^'^'^ and 

Ai(o-, f , ■!?) = V • V + (a • V)cr - eAa + ea, 
< A2{a, V, i9) = -fiAv - {fi + /u')V(V ■v) + Va - KjiVAa - Kj2^A'd + ev, 
A3{a, V, i9) = -qAi? + €'&. 

Clearly, A is closed operator. Furthermore, Lemma 2.1 implies that for each < e < eo, the range of 
A is closed. 

Proposition 2.1. There exists two positive constants 6o = ^0(71, 72, k, A*, /u', d) and eo = eo(7i) 72, k, 
fi,IJ,',a) < 1 such that if S in ^2.4^ satisfies 6 < Sq and < e < eo, then \2.T^ has a solution 
{a,v,-d) G T-L^''^'^, which satisfies 

Il(^>^>^)ll3,2,3 < C{e){\\{g,f,h)\\ + \\V{g,h)\\). (2.25) 

where C(e) > is a constant depending only on 71,72, k;,//,/^', d and e, and C(e) ^ cxo as e ^ 0. 

Proof. Firstly, for any {a,v,^) € "H^'^'^ and ((T*,'L'*,t?*) G -^00,00,00^ ^^ have from integration by 
parts that 

{A{a,v,'d),{a^,v^,'d^)) = {V ■ v + {a ■ V)a - eAa + ea,a^) + (-dAi? + ei?,T?*) 

+ {-fiAv - (/i + /u')V(V • u) + V(T - K7iVAo- - k72VAt? + ev, u*) 
= (<T, -V • u* - V • (ao-*) - eAcT* + ea* + k7iA(V • u*)) (2.26) 

+(t;, — /iAu* — (/i + /i')V(V • u*) — V(T* + ew*) 
+(t?, -dAt?* + ei?* + K72A(V • t;*)) 
Therefore, D{A*) = ?^2,3,2 ^^^^^ f^^ ^^^^y (ct^^^;^^^^) g -^2,3,2^ 

A*(cr*,t;*,'!?*) = (Ai(cr*, V*,1?*),^2('^*;^*>^*))^3('^*;^*;^*)) > 
where 



^i(cr*,t;*,'!?*) = -V • f * - V • (ao-^,) — eAcr^, + ea^ + k7iA(V • f*) 
^2('^*>'"*'''^*) = -/wAi;* — (/U + ^')V(V • V*) - V(T* + eu*, 
^3(0-*, u*,-!?*) = -aA-d^, + e-!?* + k72A(V • !>*). 



(2.27) 
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Employing the same argument as in the proof of Lemma 2.1 and using the equation: 

A(V • v^) = - ^ ^ , (Act, -eV-v^ + V-A*) 
2ii + ^' 

which follows by taking the divergence "V- " on both side of (J2.27p 9. one can get 

||A(V • v,)\\ + \\V{a„v„i},)\\i + e\\{a,,v,,^,)\\2,i,i < C,\\{Al,A*2,AlV- Al)\\ (2.28) 

Hence the closed range theorem implies the existence of solution to (|2.7p . (j2.25p follows directly from 
()2.8p . This completes the proof of Proposition 2.1. 



2.1.2 Solution to linearized problem (|2.3p and its L^ estimate 

In the following Lemma, we discuss the estimate for solution to ()2.7p independent of e. 

Lemma 2.2. Assume that (<t, f,i?) S 'H?''^'^ is a solution of the approximate problem, {2.1/\) . Then 
there exists a constant 5q = (5o(7ii725 1^, tJ-, fJ'' , o) > such that such that if 6 in i2.4\ ) satisfies 6 < Sq, 
we have the estimate 

\\V{a,v,i})h,^,5<C{\\il + \x\){g,f,h)\\ + \\V{gJ,h)h,2,3}, (2.29) 

where the constant C > depends only on 71,72, K,/U,/i' and a. 

Proof. Using the Friedrichs mollifier, we may assume that {a,v,'d) G 7^°°'°°'°°. By the same 
argument as in the proof of Lemma 2.1, we have 

1 
I|V(ct,t;,t?)|||i,2 < Ci{||/f + \\V{g,f)f + J^((V-<?, W) + (V^/, V^r;) + (V^/i, V^^)) } (2.30) 

For the third term on the right hand of (j2.30p . the Cauchy inequality and the Hardy inequality imply 

that 

1 

u=o (2.31) 

<^jV{a,v,^)\\l + C\\il + \x\){g,f,h)f 

Consequently, 

l|V(cT,t;,^)ll|i,2 < C{||(1 + \x\){g,f,h)f + \\V{g,h)f} (2.32) 

where the constant C depends only on 71, 72, k, /x, fi' and a. 

Moreover, for any multi-index a with 1 < \a\ < k — 1, applying 9" to (|2.7p . we have 

'V-d^v + {a- V)dy - eAdy + ed^a = d^g - la, 
-fiAdy -{fi + fi')V{V ■ d^v) + Vdy - K-fiVAd^a - K-i2VAdy + ed'^v = d^f, (2.33) 
-dA5°t? + ed^^ = d^h, 

where 

(3<a 
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with Co bemg the binomial coefficients corresponding to multi-indices. Notice that the third term 
on the right hand of (j2.30p can also be estimated as follows: 



i/=0 



(2.34) 



Thus, it follows from (12:30]) and (l23il) that 

||V(cT,«,T?)||i,i,2 <C{\\{gJ,h)f + \\Vig,h)f + \\{a,v,i})f} (2.35) 

Applying ([235]) to ([233]) . we have 

\\Vd^{a,v,{))\\l,^, <C{\\d^{gJ,h)f + ||Va:(5,/i)f + \\d^{a,v,^)f + \\ljj} (2.36) 

Since 

\\Ia\\l<C6^\\Va\\l^ (2.37) 

as follows from the Sobolev inequality and the assumption (|2.4p . We get from (|2.36p and (|2.37p that 

||Vl"l+V,Vl°l+2t;,Vl"l+¥f 

r -1 (2-38) 

<c{\\d-{g,f,h)f + \\Vd-{g,h)f + \\d-{a,v,^W + 6^Va\\l^} 

Combining ()2.32p and ()2.38p . we obtain (|2.29p if 5 > is small enough. This completes the proof of 
Lemma 2.2. 

Now, we are ready to show the existence of solution to the linearized problem (|2.3p by using 
([2:29]) . 



Proposition 2.2. There exists 6q = (5o(7i,72) ^)/i;M',a) > such that such that if 6 in i2.4\ ) 
satisfies 5 < 5o, then the linearized problem 112. 3\) admits a solution (17,^,1?) G 7^6,5,6 yjji^QJi satisfies 
the estimate: 

\\{cJ,v,mL^ + l|V(^,t;,^) 115,4,5 < C{\\{l + \x\){gJM\ + II V(5, /, /i) 113,2,3} (2.39) 

where the constant C > depends only on 71,72, K,//,/i' and a. 

Proof. Set 

i^=||(l + |x|)(r7,/,/i)|| + ||V(5,/,/l)||3,2,3 

From Proposition 2.1 and Lemma 2.2, it follows that for each < e < eo(eo is given in Proposition 
2.1), (1221) admits a solution (c^^^;^??^) G Ti^^^^^ which satisfies 

l|V(a%t;S#)||5,4,5<Ci^. 
The Sobolev inequality imply that 

\\{cj\v\^^)\\L.<C\\V{a\v\^^)\\<CK. 
Choosing an appropriate subsequence, there exist {a,v,'d) G L^, {ai,Vi,i}i) G 'H^'^'^ such that 

{a'^,v'^,'&'^) ^^ {a,v,'d) weakly in L^ 
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/da' dv' d-d' 
\ dxi ' dxi ' dxi 



ai,Vi,'di] weakly in Ti 



5,4,5 



as e ^ 0, then one can check easily that 

da dv &•& 



and 



OXj OXi OXj 



(a,t;,^)|| + ||V(a,7;,i9)||5,4,5<CE:. 



On the other hand, we have 

V -v' + {a- V)a' - eAa' + ea' — > V ■ v + {a ■ V)a, 
-fiAv' - (/u + /i')V(V • v') + Va' - K7iVAcr^ - k72VAt?^ + ev' 
— > -fiAv - {fi + fi')V{\/ ■v) + Va - K7iVAcr - K-f2VA'&, 
-aA'&' + ei^' — > -qAi?. 

in distribution sense. This completes the proof of Proposition 2.2. 

2.1.3 Weighted L^ estimate for solution to the linearized equation (12. 3p 

In this subsection, we give the weighted L^ estimate for the solution to (|2.3p . 



Lemma 2.3. Let {a, v, i?) G 7^^,5,6 ^g ^ gQinHon to Ii2.&]) which satisfies Ii2.&y\) . Then there exists a 
constant 5o = (5o(7i,72, K,/i, /^',a) > such that such that if 6 in \2.4^ satisfies 5 < 5o, we have for 
any 1 < / < 4 that 



Y,{\\{'^ + \x\y ( W, V+^a, V^+2^) II + II (1 + ix|)^ (V^+^w, V^^^, V^^^) II } 



u=l 



(2.40) 



<cl\\V{a,v,^)\\ + |16i|1j5||ci||j5 + ||62||,/5||c2|1jv5 + ^ 11(1 + 1x1)^^^-^/,^) 

I u=l 

+J2h + \x\rv'^{gM] 

v=\ ) 

where C > is a constant depending only on 71,72, K,/i,/i' and a. 

Proof. The proof is divided into four steps. 

Step 1. Using the Friedrichs mollifier, we may assume that {a,v,'Q) € 7^°°'°°'°°. For any multi- 
index a with |a| =1 — 1, applying 9" to (|2.3|) 9. then taking the L^ inner product with (1-|-|xJ)^'V5"(t 
on the resultant equation and summing up a, we have 



(l + |x|)'(V'a,V'+V 
< C{|(|V'+it;I, (I + \x\f^Va)\ + |(V'(7, (1 + |x|)2'-iV'+V)| 

+ |(V'A^, (l + |xj)2'vV)|+ (V'-H-(^i-V)ci + /},(1 + |xI)2'v'(t) } 



(2.41) 
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From ()2.3p 9. we also obtain 

2 



(l + |x|)'V'+V 



< C <^ 11(1 + IxlYV^af + 11(1 + IxiyV^+'^vf + 



;i + |x|)'v'-v 



+ (V'A^ + V'-H(^i • V)ci}, (1 + [x|)2'V'+2,t) 



Thus, it follows from a linear combination of ()2.4ip and (|2.42p that 

2 



lV7l+l Il2 



<Ci ||(l + |x|)'V'+^7; 



(l + |x|)^V'-V + \{\V^+^v\,{l + \x\)^^Va)\ 



+ |(VV, (1 + |x|)2'-1v'+1ct)| + |(V'A^, (1 + jx|)2'V'a)| 



+ 



(V'-i/, (1 + \x\fV^a) + |(V'A^, (1 + |x|)2'v'+M| 



+ (V'~H(^i • V)ci}, (1 + |x|)2'V'+M + (V'-H(6i • V)ci}, (1 + \x\y'V'a 



cl \\(i + \x\yv^+^vf + 



+ Il + l2+h}. 



(l + |x|)'V'-V 
The cauchy inequality implies that 

/i < ?? 11(1 + |x|)'(VV, V'+V, V'+2a)f + C^{||(1 + \x\yv^+^vf 

+ 11(1 + Ix|)'-iV'a||2 + 11(1 + |x|)'V'+¥f + 11(1 + |x|)'V'-i/f } 
For /2, notice that 

(1 + |x|)i"i+i'^i!a°6i!|afci| < c||fe||j5||c||j5 

for any multi-index a, (3 with \a\ < 1 or |/3| < 1 and |qJ, |/3| < 5. If 1 < / < 3, since 

(l + |x|yv'-H(&i-V)ci}GL2 
as follows from (|2.45p . we have 



I2<V 



_i_ r< \\h i|2 11^ i|2 



(l + |x|)'v'+V 
If / = 4, we get from integration by parts that 

h < |(VH(6i • V)ci}, (1 + \x\f'V^+^a)\ + |(V'-H(foi • V)ci}, (1 + |x|)2'-iv'+V) 

= h,! + -^2,2 

Using the Leibniz formula, we have 

h,i <Y. Y. C^c. \{{dr^h, ■ V)afci, (1 + \x\fv'-^'a) 

|a|=4/3<a, |/3|=1,2 

+ E E ^" ^(^""''^1 • v)^^i' (1 + i^i)''"'v'+V) 

|a|=4/3<a, |/3|=0,3,4 
— -'2,1 T^ -'2,1 



(2.42) 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



(2.47) 



(2.48) 
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By integration by parts, I2 1 can be estimated as follows 

4,1 < C 11(1 + |x|)2Vfei||^^ { 11(1 + |xi)2V3ci|| 11(1 + |x|)4VV| 
2 

+ J^ 11(1 + ixD'^v'^+^^ill 11(1 + ixD^vVII } 

+(the same term except for the exchang of 61 and ci) 

<r?||(l + |x|)4(V5fj,VV)||' + C^||6i||25||ci|l25 
For I21, we deduce from ()2.45p that 

III < ^ 11(1 + !^l)^vV||^ + c^ll^ill^sllciii^s. 

Combining H^LU^i - H^IM^i . we obtain 

hi < ^ 11(1 + |x|)4(vV, vMIl' + CJh\\]4ci\\%. 
Similarly, we can also get 

l2,2 < ?? 11(1 + Ix|)4(vV, V6a)f + C^ll^ill^^llcill^s. 



Thus, it follows from (ITiH]) . dHIl), (l23T]l and ([232]) that 



I2<V 



(l + |xl)'(V'+V,V'+V) 



I /^ 1 1 J, 1 1 2 II 1 1 2 
+ '-"')ll"l|lj5||Cl||j5- 



Similar to the estimate of I2, we have if 1 < / < 3, 



h<r] 



'1+ Ixh'V'cr 



_1_ r' \\h l|2 lU ||2 
"T "-^ryllfl II J5 ||l.1 II J5; 



and if Z = 4, 



n4v74,||2 



3v73,||2 



h<v 11(1 + kirvvji + c^iibiii^siiciii^s + 11(1 + kD'^vv 



Substituting ([2331), (f233D - (f235D into (JTiSll . we arrive at 

'^V7'^ V7'+l^ T7'+2, ^ 



(l + |x[)'(VV,V'+V,V'+V 

<c\\\(i + \x\yv^+^vf + 



'i + \x\yv^-^f 



+ 11(1 + |x|)'V'+ii9| 



(1 + Ix|)'-i(V'-V, VV)f + llfoill^sllcill^,} 



(2.49) 



(2.50) 
(2.51) 
(2.52) 

(2.53) 

(2.54) 
(2.55) 



(2.56) 



Step 2. For any multi- index a with \a\ = I, applying 9" to ()2.3p 9. then taking the L^ inner 
product with (1 + \x\)'^'^^^v on the resultant equation, we have from integration by parts that 



H\\{1 + |x|)'va>f + ^(V9>, 2/(1 + |xl)2'-i^a°7;) + {fi + fi')\\{l + jx|)'V • d^^v 



iO.,ll2 



+ (^ + /.')(V • 5,"z;, 2/(1 + \x\f-^j^^dy) + (Va°a, (1 + IxD^'^^t;) 
-K7i(VAa^CT, (1 + \x\f^d^v) - K-f2{VAd^^, (1 + \x\f^d^v) 
= {VAd^{-{bi ■ V)ci + /}, (1 + |x|)2'a» 



(2.57) 
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Applymg 9° to (|2.3|) i . then taking the L^ inner product with {l + \x\)'^''d^v on the resultant equation, 
we have from integration by parts that 



-{d^v, (1 + |x|)2'va-a) - {d^v, 2/(1 + \x\fi-^^dSa) 

+ (Va^((a • V)a), (1 + Ixlf'd^a) = {d^g, (1 + \x\)^'dSa) 



(2.58) 



Cancehng the term —{d^v, (1 + |x|)^'V9°(t) by adding (j2.58p to (|2.57p . and taking summation with 
respect to a, we obtain 

(1 + |x|)V'+^y ^ < C K |(V'+if, (1 + !x|)2'vM| + |(V'i;, (1 + Ix|)2'-1v'(t)| 
+ \{V^g, (1 + |x|)2'v'ct)| + |(v7, (1 + |x|)2'v't;)h 
+ |(VH(a • V)a}, (1 + |xl)2V'(j)| + |(V'{(6i • V)ci}, (1 + \x\f^V^v)\ (2-59) 

+ I (VAV'a, (1 + |2;[)2' V'u) | + | (VAV't9, (1 + \x\f^V^v) \ } 

= C{/4 + h+h + l7 + h} 

Integration by parts and the Cauchy inequahty imply that 

h<ri\\{l + \x\y{V^a,V^+\)f 

+C7^ {||(1 + \x\y-^V^vf + 11(1 + \x\yv^gf 

Similar to the estimate of /2 i, we have 



1 + b )'V' 



lV7l 



-'fY} 



h<v 



(l + |x[)'V'+^« 



+ Cr,\\\(l + \x\y~^\/'vf 



For Ij, we deduce from integration by parts and ()2.3P i that 



l7< 



(AV'cr,2/(l + |xj) 



21-1 X yil 



V'v) 



+ \{AV^a,{l + \x\y^V{V ■ v)) 



\h l|2 II ||2 \ 

Pl II jfc+i l|Ci II jft+i r 



< 21 |(V'+V, (1 + !x|)2'-i|V'f 1)1 + |(V'+2fj, (1 + [x|)2'v'{(a • V)a})\ 
+ |(V'+V,(l + |x|)2'-iV'5|)| 

< r? ||(1 + |x|)'V'+2a||^ + C^ (11(1 + \x\y-'^V^vf + 11(1 + |x|)'V'5||2) + A 



where 



A = |(V'+V, (1 + |x|)2'vH(a • V)a})| 

= J^|(v'+V,(i + !x|fa:{(a.vM) 

|a|=« 

< J^ |(V^+V, (1 + I^l)^'(a • V)9» 

|o|=« 

+ E 5^ Cf |(V'+2^, (1 + \x\fd^-^a ■ d^Va))\ = Ai + A2 

\a\=ip<a 

The Sobolev inequality and the Cauchy inequality imply that 

Ai < C\\a\\Loo ||(1 + Ix|)'V'+V|| ||(1 + |x[)'V'+V|| 
<C6\\{l + \x\y{V^+^ay+^<j)f 



(2.60) 



(2.61) 



(2.62) 



(2.63) 



(2.64) 
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H=^ /3<a,|a-/3|<ifi /3<a, |q:_/3|>M 

^E E cf||(i + N)'"-''-^^^ 

l"l=^/3<Q,|a-^|<H 

+ E E c^ 

l°l='/3<a,|a-/3|>M 

' II 2 

<C75^ (l + |x!)'(V'a,V'+V,V'+V 



l+[2;|)l°-/3|-ia"a 



i + l^l)i^i+ivi/^i+V 
(i + |x|)i'3|+ivi^i+V 



i + l^iyv'+V 



!/ = « 



Combining ()2.62p - ()2.65p . we obtain 



Ir < 7? 11(1 + \x\yV^+^af + Cri (||(1 + \x\y-^Vvf + ||(1 + \x\)^V^gf 



+C5Y,\{'^ + kl)'( vV, V'+V, V+V) 



(2.65) 



(2.66) 



Finally, similar to the estimate of A and Ij^ respectively, we have 

I 



h<C6Y, (1 + I^I)'(V'ct, v'+V, V'+V) 



!/=! 



(2.67) 



and 



Is < 77 11(1 + Ix|)'V'+2^||^ + Cr, (\\(l + \x\y~^V^vf + 11(1 + IxlYV^qf 



+C5^\\{1 + |x!)'(V'cT, V'+1(T, V'+^T?) 



1^=1 



(2.68) 



Substituting i^Wh . i^Mh . (IXUHD - dTUgD into (ITMIl . by the smallness of 7/, we arrive at 



iv7«^l|2 



(1 + |x|)'V'+^ < C |7? 11(1 + \x\yVaf + {r] + 6) ||(1 + |x|)'+^V'+^(a,i?)|| 

+C7^f 11(1 + IxD'V'sll^ + 11(1 + kD'-^V't;!!^ + (1 + k|)'V'-i/ 



(2.69) 



+^5^||(l + |x|)'^(VV,V'^+V)||'} 



Step 3. For any multi-index a with |a| = I, applying d" to ()2.3P '^. then taking the L^ inner 
product with (1 + |x|)^'(9"7? on the resultant equation, integrating by parts and summing up a, we 
can get 

Q 11(1 + \x\yV^+H\f < C |(V'+^7?, (1 + |xj)2'-lV'7?|)| + C {V^{-ib2 ■ V)C2 + h}, (1 + [x|)2'v'7?) 

= l9 + ho- 



For /g, the Cauchy inequality imply that 



a 



(1 + |x|)'V'+^7? 



+ c 



'i + \x\y-'v''d 



(2.70) 
(2.71) 
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Shnilar to the estimate of Jy, /lo can be estimated as follows 



a 
ho<-r 



(l + |x|)'V'+^^ 



+C7<^ (l + |x!)'-iV'^ +||62r75|lc2 



\N' 



+ 



Putting (f2TT1) and (l2:72]l into l^7fU\i gives 



IV71+1. 



{i + \x\yv'+'^ 



<c 



l-l^L 



'i + \x\y-'Vi} 



+ Il^2||j5|lc2||jv5 + 



On the other hand, we also get from ([231)3 that 



(1 + |x|)'V'+'t? 



<c 



<C[ ||62||j5||c2||^5 + 



l + lxh'V'-^/i 



;i + [x|)V-i/i 



;i + |xj)V'/i 



(2.72) 



(l + |x|)V/i 

Consequently, we deduce from (|2.73p and (j2.74p that 

||(l + |x|)'(V'+i^,V'+2^)||^ 

<c|||(l + jx|)'-iV'^||' + |162|1^5||c2||^5+ {l + \x\)\V^-%V^h) ' 



(2.73) 



(2.74) 



(2.75) 



Step 4. Now, we begin to prove (|2.4U|) by using the estimates in the above three steps. We use 
the method of induction. First, for the case of / = 1, we derive froin (|2.43p . ()2.44p . (|2.46p and ()2.54p 
that 

- 11(1 + \x\)V^'d\\^ 



||(l+Ix|)(Va,VV,VV)||^ <C\ ||(l + |a;|)V2?;||^+ (l + lxj)/ 

+ \\V(y\\^ + \\hifjrAci\\%Y 
which, together with (|2.69p and (|2.75p with / = 1 gives 

\\{l + \x\){Va,V^a,V^a)f + \\{l + \x\)V^vf + \\{l + \x\){V'^^,V^M^ 



(2.76) 



< C \\V{<J,V,m + ll^llljsllcillJB + I|fe2|l55||c2||^5 + (1+ Ix|)(/,/l,V5,V/l) 



(2.77) 



by the smallness of rj and 5. Thus, we assume for / > 2 that 

;i + |x|)'-i(V'-V,V'a,V'+V) V {l + \x\y-Wv V 
<c[\\V{a,v,d)f + llftill^sllcill^s + ||62||^5||C2||^5 

+ ^||(l + |x|)'^(V--i/,V"/i,V'^5,V'^-i/i)| }. 



{l + \x\Y-'{V^,V+'{}) 



«+i. 



(2.78) 



v=l 



Furthermore, the linear combination [Mi x ([ZBO]) + ([236]) ] + M2 x ([2775]) for Mi > and M2 > 
large enough in turn gives 



l + |x[)' VV, V'+ V, V'+^CT 



+ 



(1 + Ix|)' V'+^77,V'+^T?,V'+^^ 



< C 11(1 + |x|)'-i(V'-V, V'a,V't;, V'^)|P + ||6i||25|lci|l25 + I162fj5||c2||^5 



(2.79) 



+ 



i-i 

(1 + ix|)'v'-i(/, h)\ + ||(i + ix|)'v'(<7, h)\ + 6Y,W + \x\r (VV, V^+V) II' } 



u=l 
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provided that r/ and 5 are small enough. Combining (j2.78p with (|2.79p . if 5 > is small enough, we 
can get (j2.40p . This completes the proof of Lemma 2.3. 

Combining Proposition 2.2 and Lemma 2.3, we have the following theorem. 



Theorem 2.1. There exists 6o = (5o(7i,72, k, /i, /u', a) > such that such that if 6 in \2.4\) satisfies 
5 < 5o, then ^2.'J\) admits a solution {a,v,'i}) G 7^6,5,6 yjjiiQJi satisfies the estimate: 



(2.80) 



||(a,z;,^)||i6+^ 11(1 + |x|)'^(VV,V-+V,V'^+V)||+^ 11(1 + |x|)'^-iV'^^|| 

iy=l u=l 

5 

+ ^||(l + |x|)'^-i(V'^i?,V'^+i^)|| 

v=l 

< C'{ll&illj5||ci||j5 + ||62|1j5||c2|U5 + 11(1 + \x\){g,h)\\ 

4 3 

+ Y^\^(l + \x\YV-{g,h)\+Y,\{^+\Ar^'^''Cf,h)\] 
where C > is a constant depending only on 71,72, k, //,//' and a. 

2.2 Proof of Theorem 1.1 

In this subsection, we shall construct a solution to ()1.6p by the contraction mapping principle in 
Ae' ' . To this end, we consider the following iteration system 

f Op 

V-i; + — (t;- V)cr = g, 

P . . 

-fiAv - (/u + /i')V(V ■v) + Va- KjiVAa - K72VA1? = -p{v ■ V)v + /, (2-81) 



^ -qAi? = -fji{v ■ V)i? + h, 



where 



g = — ^(u- V)t?H ^, 

P P 

f = -{p - p){v ■ V)v + Hip (Va ■ V^pp + Vpp ■ V^a + VppAa) + k {ppp - ppp) VAa 



+KP [y^ ■ V'^pe + Vpe ■ V^T? + VpeA^j + k {ppg - ppg) VAi9 + pF - vG, .^ ^^. 

h = -{m- m){v ■ V)i? + miv- V)a + ^{v) -fi^G + H + $(p, v) + —G- GyOG, 

^i = /5C7v-M n. = '-^. m=~^. = 9 + ^. 
ppp p ppp 

Here, (cr, {;,■!?) G A^' ' is given, and pp = pp{P + a,6 + '&),fii = r]i{P,6), etc. 

2.2.1 Introduction of solution map T for (12. 1|) 

Firstly, we apply Theorem 2.1 to (j2.8ip to get the weighted L^ estimate. Let 

PP _i _ ~ 

a = — —V, 6i = ci=p2-u, 1)2 = fill), C2 = t?, (2.83) 
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and g,f,h m Theorem 2.1 be defined as in (j2.82p . We choose e > sufficiently small such that 
^ < p < 2p, as follows from the sobolev inequality. Assume that the assumptions of Theorem 2.1 
hold and denote 

3 

Ko = ^ 11(1 + Ixir+^ViG, F, H)\\ + 11(1 + \x\)^V\G, H)\\<oo (2.84) 

then we can check (|2.5p and ()2.6p easily and additionally we have 

4 3 

11(1 + \x\){g, h)\\ + ^(1 + IxirV'ig, h)\\ + ^(1 + Ixir+'V'if, h)\\<C (e^ + K^) . 

for some constant C = C{p,9, fi, p' , k) > 0. Applying Theorem 2.1 to (|2.8ip . we have the following 
lemma. 



Lemma 2.4. Let {G,F,H) G H ''^' satisfy 1^2. S4 )- Then there exists a constant eo > such that if 



e < eo, 112.81]) with {a,v,i!}) € A^' ' admits a solution {a,v,'d) £ 7^^,5,6 y^j^j^QJ^ satisfies the estimate: 

4 5 

II (a, V, mL^ + E 11(1 + i^D' (^''^' ^'^'^' ^'^'^) II + E 11(1 + 1^1)'^' v^HI 

5 ^^=1 ^^=1 (2.85) 

+ ^ II (1 + \x\Y-^ (V^^, V'^+^^) II <C{e^ + Ko) 

v=l 

where the constant C depends only on p, 9, p, p', k and a. 

Based on Lemma 2.4, we can define the solution map T : A^' ' — >■ T^^'^'^ by (a, v, i9) = T{d, v, i}). 
Since the contraction mapping principle will be applied to prove Theorem 1.1, we have to show that 
T{a,v,'d) = (iT, u,t?) S Ae' ' . To this end, we first cite the following lemma which will play an 
important role when we estimate the solution by the L°° norm. 

Lemma 2.5. (16]) Let E[x) he a scalar function satisfying 

\d-E{x)\<-^, |a|=0,l,2. 

(i) If (j){x) is a smooth scalar function of the form (j) = V ■ (pi + (j)2 satisfying 

Li{^) = \\{1 + \x\f<P\\Loo + 11(1 + kl)^0i|Loo + \\<P2\\l^ < oo, 
then we have for any multi-index a with \a\ = 0, 1 

|a,"(i?*<^)(x)|<-^Li(0). 
(ii) If (j){x) is a smooth scalar function of the form (j) = (j)i(j)2 satisfying 

L2{4>) = 11(1 + \x\f4>\\Lo. + 11(1 + k|)'(V(Al)(A2|Loo + 11(1 + \x\fM^<t>2)\\^ < OO, 

then we have for any multi-index a with \a\ = 1,2 

\d^{E*c^){x)\<-^L2{P). 
\xy"\ 

Here Ca denotes a constant depending only on a. 
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With the aid of the Helmholtz decomposition and the Fourier transform, the solution of (j2.8ip 
can be formulated as follows, cf. |6|. 



where 



V = w + Vp, a - K-yiAa = $ + K'y2A-& + (2/i + fj,')Ap, i^ = E*Q, 



i=l 

p{x) = Eo*Rix), 
I ^ = £^0 * (V • /). 



(2.86) 



(2.87) 



and 



Eij{x 



1 / Oij XiXj 



STTfJ, \\x\ \x\'^ 

fi = -p{v ■ V)Vi + fi, 
R{x) = -^{v-V)a + g, 

e = -{fii(v-v)^-h}. 

a 



En 



47r rr 



(2.88) 



Now, we shall estimate the L°° norm of the solution to (j2.81|) by using Lemma 2.5. 
Lemma 2.6. Let (G, F, H) G //^'^.^ satisfy the following estimate: 

K ^ \\{G,F,H)\\c + 11(1 + \x\fv\G,H)\\ < oo 



(2.89) 



// {a,v,i!}) € 7^^'^'^ is a solution to i2.81]) with (a-,{i,i9) G A^' ' and satisfies i2.85]) . then (a, f,i9) 
satisfies the estimate: 



j; 11(1 + \x\fV^a\\^^ + E 11(1 + kir'V^(«,^)| 
where the constant G > depends only on p, 9, p, p' , n and a. 



L°° 



(2.90) 



,4,5,5 



Proof. First, we deduce an estimate on /. Since (a,?),'!?) € A^' ' , there exits Vi = (V]*)i<j<3 and 
V2 such that V ■ V = V ■ Vi + V2, and 



{l + \x\fVi 



+ 



{l + \x\)-'V, 



<e 



Thus we have 

fi = -p{v ■ \7)vi + Kp[Va ■ V'^pp + Vpp ■ V^o- + VppAa]i + k {ppe - ppe) AOx, 

+Kp[V{} ■ V^pe + Vpe ■ V^ + VpeA^i + k {ppp - ppp) Aa^^ + pFi - iiG 

= V • (^-pviv + pviVi + pFi,i^ + {-p(Vl • V)vi - Vi{Vi ■ V)p + pviV2 + pF2,i 

+Kp[Vi} ■ S/^pg + S/pe ■ V^i^ + S/peA^]i + Kp[S/a ■ V^pp + Vpp ■ S/^a + S/ppAa]i 
+Vi{v ■ V)p - Vp ■ Fi^i + K {ppp - ppp) Aa^^ - ViG + k {ppe - ppe) A^^^^ 
= V • Ai + /, 2. 
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Here [•••]{ denotes the i — th component of the vector [•••]. 

(j2.85p . using the Sobolev ine 

1(1 + \x\ffiL^ + 11(1 + I^DViJU + Whih^ < C{e^ + Ki 



By (cjjt;,'!?) € A^' ' and (j2.85p . using the Sobolev inequahty and mean value theorem, we have 



and 

11(1 + IxDM^ILoo + 11(1 + \x\)^M\l^ < C{^ + ifi) 

where 

Ki = 11(1 + |x[)3(F,G, VF,VG)||^^ + 11(1 + IxD^Fill^^ + | 

Hence, by (i) and (ii) of Lemma 2.5, we obtain 



[x|jt(;j|, \x 



|$|, |V$|, \Vwj\, \V'^Wj\) < C{e^ + Ki). (2.91) 



As for Vp, V^p, V^p, due to [7], 

|x||Vp|, |x|2 (IV^pI, |V=^p|) < C{e^ + Kq + K2). (2.92) 

where 

K2 = 11(1 + \x\)^G\\^^. + 11(1 + |x|)3(VG, V^Oll^^ 

Combining dSSSDi, (|2:9n) and ([2:92]) yields 

[x|[t;[, \x\^{\Vv\,\V^v\) <C{e^ + Ko + Ki + K2). (2.93) 

Next, we turn to estimate "d. To this end, we rewrite Q as 

e = -\miv-V)^-fl2{v- V)a - ^{v) + f,^G-H- $(p, t)) - — G + Cv(t? + e)G 
= ^V • <^ im^ - md){v - Vi) + mGi -Hi- -Gi + Cy^Gi 

H — V ■ \(Vi ■ V)(t?it9 -7?2<7) - ivi'^ -'n2a)V2 - Vfuv'd + Vfi2va - ^(u) - $(p,i)) 
a y. 

v'^ ~ ~ - 

-Vr/3 • Gi + 7?3G2 + H2 + V-vv-Gi- —G2 + CyV^ • Gi - GvG2(t? + 0) 

= v-ei + G2. 

and 

3 1 1 r „--,2 



1 1 r ~2 ^ 

e = Y,{-^mvi)a^., + ^ r/i({) • V)^ - ^(i)) + ^G - F - <D(p,^) - ^G + Gv(^ + ^)G 

3 

Since (a,!','!?) S A^' ' , it follows from (j2.85p and the Sobolev inequality that 

11(1 + \x\fe\\^^ + 11(1 + [xD^eill^^ + 116211^1 < C{e^ + Ks), 



[1 + [xD^veall^^ + 11(1 + IxlfOsW^^ < G(e2 + K3) 
1^^ + 11(1 + \x\fiVe\)Ql\\^^ + 11(1 + |x!)3Gi(Vu,,||^„ 



[1 + \x\fe\ei\\,^ + 11(1 + |x|)3(VGi)ey _ + ||(i + \x\fe\{vei)L^ < Ge\ 
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where 

Ks = \\{l + \x\f{G,H,VG,VH)\\^^ + \\{l + \x\f{Gi,H,)\\^^ + 11(^2,^2)11^1. 
Thus, it follows from (i) and (ii) of Lemma 2.5 that 

|xp|, |x|2(lVi?|,|V^i?|) <C7(e2 + i^3). (2.94) 

Finally, for the estimate of ■!?, taking the Fourier transform on both side of ()2.86p 9. we have 

1 f f r -t-^ -3 ^ ^ ^ ^ 



f- / f-OO |j,|2 \ 

/ ( / e"*"^^^t~2 dt J ($ + K72A?? + (2/i + ^')^P)(^ - y) dy 



(47rK7i)3/2 7^3 VJo 
where L{^,p, -i?) = <I> + K72A19 + (2;U + /x')Ap and we have used the fact that (cf. |18j ) 

(1 \ 1 Z"*^^ lap 
— - ) = — - / e"*"3^^t-2 dt (2.96) 
l + '^TilClV (2AC7i)3/2yo ^ ^ 

Note that the right hand of (|2.96p is the so-called Bessel potential. We deduce from ()2.9ip . ()2.92p 
and p:MD that 

N k(^ < 7 — ,,,0 / / e 4K7i*t-2dt ^ ^dy 

' ' ' ^ ^'- (47rAC7i)=^/2 ^3 vyo y |a:-yP 

X sup {Ix -yp |$ + K72Ai9 + (2/x + //')Ap| (x -y)} (2-97) 

(2:-j/)eK3 

< C7(e2 + K)-A 
where K is defined by (j2.89p and 

^= / / e 4^71*^-2 dt , ' , dy 

Jr3 Vjo / k - yr 

00 3 / /■ ixP _/__mL \ 

e"**~2 / ' ' e *-^i* dy di 
VJ]R3 \x - vY J 



roo / r 11-^ \ /"oo / r I |2 i |2 \ 

<2/ e'H~^ e'^'dy]dt + 2 e'H~^ ( , '^' ,^ e~*~^dy | dt 
Jo VJr3 y 7o \Jr3 \x - yV J 

f°° f /"°° -^ I r l7;|2 lal^ \ 

< 16(k7i)3/2 / e-*dt / e"l^l'dz + 2/ e-H'^ [ '^' ,„ e~^^^^dy dt 



00 ^ I r U,|2 

"2' 



/b(x,i) |a:-yi y (2.98) 



+2 / e"*t"t / , '^' ,^ e ^^^^i' dy ] di 



WlR3\B(^_i) |x- y| 



00 , / /• 1 



<C + 8k7i/ e^*t"2 I / ^ — ^jyUt 



yjB(x, 1) F — y 

t^-l / / L.|2, ^ 



+2 / e~'t"2 / |y|^e ''«^i* dy dt 

Jo \Jr3\B{x, 1) J 

<C 
Here -B(x, 1) denotes the unit ball in M'^. Consequently, it follows from (j2.97p and (|2.98p that 

\x\'^\a{x)\<C{e'^ + K) (2.99) 
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Differentiating the equation (|2.95p and notice that 

II IxplV^^m^^ < C ||V(Ix|2|v3^i)||^ < C{e^ + Kq), 

by using the same argument as above, we can also obtain 

|x|2|Vct(x)| < C{€^ + K) (2.100) 

Next, we consider the case of \x\ < 1. The Sobolev inequahty and ()2.85p imply that 

\\ia,v,^)\\Loo<C\\Via,v,^)t<C{e^ + Ko) 

2.101 
\\V''ia,v,^)\\Loo<C\\V''+\a,v,^)\\i<Cie^ + Ko), i, = l,2. 

(IXUnjl thus follows from (gMD, (1231, ^M, (I^^TOUjl and (IXTUni . This completes the proof of 
Lemma 2.6. 

In the following Proposition, we show that the solution (o", v,"!?) € A^' ' . 

Proposition 2.3. There exits cq > such that for any sufficiently small constant e > 0, if 
{G, F, H) G ?^4'3'4 satisfies 

K + II (1 + |x|)^ ^ILi — '^oe {K is defined in Lemma 2.6) 

then \2. 81\) with {a^v^i)) G A^' ' admits a solution {a,v,'d) = T{a,v,'&) G A^' ' . 

Proof. By Lemmas 2.4 and 2.6, it follows that ()2.8ip has a solution {a, v, i?) G 7^^'^'^ , which satisfies 

II {a, V, ^) 11^4,5,5 < C{e^ + K)< C{e^ + cqc), 

where the constant C > depends only on p, 9, fx, jj,', k and a. Thus if we take cq < -^ and e > is 
small enough, it follows that (o", f,t?) G A^' ' . Finally, we define Vi and V2 by 

V\ = va, V2 = V • — V a v ■ Vu H 

P \ P J P P 

Then it follows from (|2.8ip i that 

y .v = V ■V1 + V2. 

Moreover, by ((T,-u,t?) G A^' ' , (j2.85p and ()2.90p . we have form the Sobolev inequality that 

11(1 + IxDVill^^ + 11(1 + |xj)-V2||^, <C{e' + K+\\{l + |x|)-iG||^i} 

< C(e2 + coe) < e 

if Co < 2C ^'^^ e > is sufficiently small. This completes the proof of Proposition 2.3. 

2.2.2 Contraction of the solution map T 

In this subsection, we shall show that the solution map T for ()2.8ip is contractive. Suppose that 
(aJ,{i^7?J) G Ap'^ and (a^ , u^' , i?^' ) = T{a^,v^,^^) for j = 1,2, then we deduce from ^im\i that 

P 

-fiA{v^ - v^) - (^ + //')V(V • [v^ - v'^)) + \/{a^ - a^) - K-fiVA{a^ - a^) (2 102) 

-K72VA(t?l - 1?2) = _p2(^l _ ^2) . y~l _ -2(-2 . y)(^l _ ~2) ^ f^ 



^ -aA{^^ - t?2) = -f^l[{v^ - v'^) ■ V¥ + {v^ ■ V){¥ - T?2)] + h. 
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where 



P^ P^ 



7d 



Va2 - ( ^^1 • V^^ - ^{^2 . ^^2 j ^ / ^ 



1 1 



G, 



/ = -(p^ - p^){v^ ■ V)^i - (pi - p2)F - (-^i - v'^)G + [kp^ (Va^ • V^~p\.+ 

+Vp\, ■ V^^i + Vp],Ad^) - Kp^ {Va^ ■ V^pI + Vpl ■ V^d^ + Vp^Aa^) } 

+ k|p1(Vt?1 • V^pI + Vp] ■ V^T?! + Vp]A^^) - p2(V^2 . y2~2 ^ y ~2 . y2^2 



+Vp2A^2) I + ^ I ( ~1 ~1 _ ^^^) vA^i - (p2p2 _ ^p,) VA^2| 
+K { {p^p\ - ppp) VA^i - {-p^pl - ppp) VA,^2} , 
h = -{f^\ - fjl){v^ ■ V)^i + {fjl - r?2)(i}i . V)a^ + 7?2((t}i . V)ai - {v^ ■ V)a^) 

+$(i}l) _ ^(^2) ^ $(pl, ^)l) _ q>{p'^,v'^) + ^{V^ + i)2)(i}l - ^2)(- 



(2.103) 



m = p'Cs, 



~eKp^e? ~j_^'pI ~3 ^'ip's'' 



p^pp 



V2 



m 



p'p'p 



9^ = e + ^^, j = i,2. 



Since 



(1 + \x\){g, h)\\ + Y, 11(1 + l^lY'^'ia, h)\\+Yl 11(1 + 1^1)'^' V'(/' ^) 



u=l 



< C{e + K) \\{a^ - a^v^ - v^,^^ - t?2) 



u=0 



A4,5,5 



(2.104) 



as follows from the Sobolev inequality for K defined in ()2.89p . Applying Theorem 2.1 to (|2.102p . we 
obtain 

5 

\\{a^ - a'',v^ - v^,^^ - ^^)\\Le + Y,\\{1 + Ixiy-^Viv' - v^)\\ 

4 

- j; II (1 + Ixl)'' {\l\a^ - a^),V''+Ha^ - a^),V^+\a^ - a^)) \\ 
=1 

5 

+ ^ II (1 + \x\Y-^ (v^(t?i - ^'^),v^+^{^^ - ??2)) II 



v=l 



(2.105) 



u=l 



<C{e + K) 



((ji -a2,i}i- 7)2^^1 _ ^2) 



A4>5,5 



Similarly, by the same argument as in the proof Lemma 2.6, we can get 
1 1 

^ 11(1 + \x\fV''{a^ - a^)\\^^ + E 11(1 + \Ar^^^''{v^ - v\ ^' - ^')| 

=0 iy=0 

+ \\{l + \x\fv\v^-v^i9^-^^)\\^^ 
<C{e + K) (a^-a^v^-i^,^^-^^) 



i/=0 



(2.106) 



+Ce{||(l + |x|)3(y/-y^2 



L°° 



+ 



A4,5,5 



(i + \x\)-Hvi - vi)\\J 



where y/, V^, j = 1, 2 are functions satisfying 



V -v^ = V -V^ + V^, {l + \x\fV^ + {l + \x\)-^vi 



< e 



(2.107) 
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Moreover, if we define V-^ ,V2 , j = 1, 2 as 

VI = -^v^a^, Vi = V • (^v^] a^ - ^v^ ■ V'&^ + - 
p^ \ p^ I p^ p^ 



(2.f08) 



then we get from (I2.102p i that 

V • (t;i - r;2) = V • {Vl - V^^) + Vi - V^ 

and 

11(1 + \x\f{Vl - V,')\\^^ + 11(1 + \x\)-\Vi - Vi)\\^, 

< C (e + 11(1 + \x\)-^G\\^,) {a^ - a^,v^ - v\^^ - ^^) 
Combining (|2.105p - (|2.109p . we obtain 

+ 11(1 + \x\)Hvl - v,^)\\^^ + 11(1 + \x\y\vi - vi)\\^, 

<C{e + K) {a^-a^,v^-v'^,¥-^'^) 



(2.109) 



A4,5,5 



(2.110) 



A4,5,5 



+Ce 



(1 + |x|)3(y/ - v,^^^ + ||(i + \x\yHvi - vi)\\^^} 



Therefore, we have the fohowing Proposition. 

Proposition 2.4. There exits a constant cq > such that for any sufficiently small constant e > 0, 
if (G, F, H) e ?^^'3'4 satisfies 

K +'\{1 + \x\)~ GLj <coe {K is defined in Lemma 2.6), 
the for (o"-' , U-' ,'!?•' ) G A^' ' and {a^,v^,'d^) = T{a^ ,v^ ,i}^),j=l,2, we have the following estimates 



\{a^-a^,v^-v\^^-i}^) 



|A4.5,5 



+ 11(1 + |x|)3(y/ - y,2)|| + 11(1 + \x\)-^vi - vi 



2^llLi 



< 



H 



+ 



(^l_^2^~l_-2^^1_^2) 

(1 + \x\fiVl - V^ 



(2.111) 



A4,5.5 



L°= 



+ 



{i + \x\r\v^-vi)\^^] 



where {¥(,¥1), j = 1,2 satisfy [2J07\ ) and (V/,K/), J = 1,2 are defined by ^KTM) . 



Hence, by Propositions 2.3 and 2.4, the contraction mapping principle imphes the existence and 
uniqueness of solution to (|1.6p . This completes the proof of Theorem 1.1. 

3 Non-stationary problem 

In this section, we consider the stability of the Stationary solution with respect to the initial distur- 
bance (po) ^Oj ''^o) • Fix p, 6 to be positive constants and let F, G, H be small in the sense of Theorem 
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1.1. We denote the corresponding stationary solution obtained in Theorem 1.1 by {P* ,v* ,9*) , and 
set p* = p + a* = p{P*, 0*). Then by direct calculations, we have the following estimate for a*: 



5 1 

L° 



k*lljv.^ = E 11(1 + kl)'~HW*, V"+V*)|| + J] 11(1 + |x|)"+i W*||^^ + 11(1 + \x\fV^a* 

u=l u=0 

<Ce, 



where the constant C > is depending only on p and 6. Thus, we have 

\\{a*,v*,0*)\\^,,s,. ^ \\a*ys + \\v*\\j, + liriUs < (C + l)e. 

For simplicity, we assume in this section that \\{a*,v*,9*)\\jr5,5,5 < e for e sufficiently small. Define 
the new variables 

a{t,x) = p{t,x) — p*, 'w{t,x) = v{t,x) — V*, 'd{t,x) = 6{t,x) — 6*, 
then the initial value problem (|1.3|) . (|1.4|) is reformulated as 
' at{t) + V • {{p* + a{t))w{t)} = -V • {v*a{t)) , 

< wit)t - jr [pAwit) + ip + p')ViV ■ wit))] + A{t)Vait) - KVAa{t) + B(t)Vi?(t) = f{t), (3.1) 
_ T?t(t) - aD*M{t) + E{t)V ■ w{t) = h{t), 
with initial date 

(a, w, i}){t, x)\t=o = (do, wo,^o){x) = {a-a*,v- v\e - r )(0, x). (3.2) 

where 

/(t) = -{v* ■ V)w{t) - {w{t) ■ V){v* + w{t)) - {A{t) - A*) Vp* - {B{t) - B*) Vr 

/i(t) = -{v* ■ V)^{t) - {w{t) ■ V) {9* + T?(t)) + a {D{t) - D*) A{9* + ^{t)) + (D(t) - D*) H 
+ {D{t) - D*) i^{v*) + <^ip*,v*)) + D{t) [^{v){t) + $(p, v)it) - ^{v*) - <l>ip*,v*)] 

+^ [D{t)v^{t) - D*v*^] G-Cy [D{t)9{t) - D*9*] G - {E{t) - E*)\/ ■ v*, 



P p Gyp Gyp 

with A* = A{p*,9*),A{t) = A[p* + a{t),9* + ^{t)) , etc. Moreover, we set Ai{t),Bi{t),Ei{t),i = 1,2 
and Di{t) to be functions satisfying: A{t) - A* = Ai{t)a{t) + A2{t)'d{t),B{t) - B* = Bi{t)a{t) + 
B2{t)'&{t), E{t) -E* = Ei{t)a{t) + E2{t)'d{t) and D{t) - D* = Di{t)a{t), respectively. 

The aim of this section is to prove Theorem 1.2. The proof consists of the following two steps. 
The first one is the local existence result: 

Proposition 3.1. Suppose that {a,'w,'&){0) G Ti'^'^'^. Then there exists a constant to > such that 
the initial value problem ll3.1\ )- [3^) admits a unique solution {a,w,'d){0) E C(0, toj ^^''^'^)- Moreover, 
(o", w, "!?) (t) satisfies 



{a,w,^){t)\\l,^^<2\\{a,w,^){0)\\l^^. 



for any t G [0, to] • 
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And the other is an a priori estimate: 

Proposition 3.2. Let [a,w,'i!}){0) G C{0,ti;'H^'^'^) be a solution to the initial value problem 113.1]) - 
13. 0^) for some positive constant ti. Then there exists a constant eo > such that if e < cq and 
supo<i<ti ||(o-,'«^,^)(i)||4,3,3, \\{a,w,T^)\\jr5,5,5 < e, then it holds 



\\{a,w,m)\\U3+ / ||(Va,Vu;,V^)(s)||l_3,3d5<C||(a,u;,^)(0)||^_3_3 (3.3) 

for any t G [0, ti], where the constant C > is depending only on p,6, fi, fi',K, and a. 

For the proof of the local existence, we can apply the H. Hattori-D. Li [12j method directly. So 
we shall devote the following sections to the proof of Proposition 3.2. 

Before proving the a priori estimate (|3.3p . let us introduce the absolute constant e > such that 
Cqc = l/4min{p, 6}, where Cq is the constant which appears in the inequality || • H^^oo < Co|| • ||2- In 
the following lemmas and their proofs, the small constant e is at least taken in such a way that 



so that 



||(a,i(;,i9)(t)||4,3,3, II (a, u;,^) 11^5.5.5 < e < e 
{p,e) e S{p,9) = |(p,0)|| <p<^,^-<p<^ 



3.1 Some estimates for f{t),h(t) and their derivatives 

Lemma 3.1. Let {a,w,'d){t) and {a* ,w* ,i}*) be satisfying 

\\{a,w,^){t)\U,3,3, II (a, t/;,t?) 11^5,5,5 < e. 

Then for a multi-index a with < |a| < 3, we have 
(i) If we write d^f{t) of the form 



d^f{t) 



a(t) 



[pAd^w{t) + {p + p')V (V • d'^wit))] + F^{t), 



p*{p* + a{t)) 

then Fa{t) satisfies the estimate: 

' \Vv*\\w{t)\ + {\v*\ + \w{t)\)\Vw{t)\ + (|V(7*| + \Ve*\ + \V^v*\) \a{t)\ 
+ (|Va*| + jVri) |^(t)| + {\w{t)\ + \a{t)\) |G1, i/a = 0, 

|q| + 1 |a| + l 

|Vl"l+V| |u;(t)| + Y^ \^''w{t)\ + Yl (l^''^*l + l^''^*! + \^'"^^v*\) \a{t)\ 
Fait)<C{ , , ^ -=i , -=i (3.4) 

|a|+l \a\ 

I/ = l u = l 

H 
+iMt)\ + |a(t)|) Y IV'^GI + 4"l(t), if \a\ = 1, 2,3. 

Here, Rp{t) = 0,k = 1,2 and Rp{t) satisfies 



\R 



Um^-i<Ce\\iV'a,V'w){t)\\^^^ 



(3.5) 
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(ii) If we write d^h{t) of the form 

d^h{t) = aDi{t)a{t)Ad^'&{t) + H^{t), 

then Ha{t) satisfies the estimate: 

' (iV^ri + |Vw*| + |Vw*||V2(T(t)| + |Vii;(t)||V2CT*|) \o(t)\ + \v*\\Vd{t)\ 
+ (|V??(t)| + |Vr I) \w{t)\ + (|Vct*| + |VCT(t)|) |VCT(t)| + |??(t)||Vt;*| 
+ (|V(T(t)| + |W(t)| + |Va*| + |VV*| + |Vu;(t)| + |Vv*|) \Vw{t)\ 
+ |a(i)||i?| + |(a,«;,i?)(t)||G|, i/Q = 0, 

|«|+1 |a|+2 



i?a(t) < C < 



|a|+l 



|a| + l 



i/=l 



+ ^ (|V"«;(t)| + |V"^(t)|) + |^(t)| Y. IV^r;*! + [w;(t)|| Vl^l+^r | 



i/=i 



|a| 



v=l 



+ |(a,«;,T?)(i)|^|V^G|+i?gl(t), i/|a| = 1,2,3. 

iJere, -R}^(t) = and Rjjit), R%{t) satisfies 

\\R%it)\\^s<Ce\\vMt)\\, WRlm^i <Ce\\iVa,Vw)m3,2 

Proof. By the Leibniz formula and the Sobolev embedding: H^ M- L°°, we can check ([37 
with 

ri?^(i) = o, if k = 1,2, 

\ R^t) = \V^wit)\ |V3cj(t)| + \V^w{t)f + I VV*| \Va{t)\, 
and 



(3.6) 



(3.7) 
(1331) 



(3.8) 



< i?3^(t) = (|v2^(i)| + |v2u;(t)|) |v2u;(t)| + (I W(t)| + |vV(t)| + \vMt)\) |v2u;(t)| (3.9) 
^ + I W(t)| |v3tt;(t)| + (I VV(t)| + I W*(t)|) |Vw(t)I + (|Va(i)I + | W(f)|) |VV| 

For the proof of (j3.5p and (j3.7p . we only give here the estimate of the most difficult term R^. The 
others can be dealt with similarly. Using the Gagliard-Nirenberg inequality, we have 

WRirn^^ < c{ ||v2^(t)||^s \\^Mt)\\ + W'^Mt^Le ||v2t«(t)|| + ||v-V(t)|L6 llvMOll 

+ ||V2u;(t)||^6 ||V3^(t)|| + ||V2u;(i)||^6 || W(t)|| + \\\/^ait)\\^, ||V3u;(t)|| 
+ ||(VS*,VV)||^6||Vu;(t)!| + ||(V(j,W)(t)||||VV||^e} 
<Ce\\{S/a,S/w)m3,2 
which is the desired estimate (jS.ZP o. This completes the proof of Lemma 3.1. 



30 



Zhengzheng Chen and Huijiang Zhao 



3.2 Estimates for Vw{t),V'd{t) and their derivatives up to V^w{t),'V^'&{t) 

Lemma 3.2. Let {a,w,'d){t) E C{0,ti;'H^'^'^) be a solution to the Cauchy problem l[3J ]) - ![3^) for 
some positive constant ti. Then there exists four constants eo, Aq > and di,d2 > such that if 
< e < eo and \\{cr,w,'&){t)\\4^3^3, \\{a,w,'d)\\jr5,5,5 < e, then it holds 



d_ 
di 



+ (i(t)Vcj(t), Vcj(t)) + {A{t)wit),wit)) + {B{t)^{t),^t)) 



+di\\Vw{t)W^ + (i2||V^(t)||2 < Ce\\Va{t)\\l 
and for 1 < k < 3 and any A with < A < Aq, 

2 



(3.10) 



d 
di 



vV(i) + (i(t)V'=+V(t), V'=+V(i)) + (i(i)V*=u'(t), V^u'(t)) 
+(S(t)V'=i9(t),V^'i9(t))) +di||V'=+iu;(t)||2 + d2||V^'+ii9(t)||2 
<C{e + X)\\V{a,w,m)\\l+i,k-i,k-i + CX~HvHw,m)r 
where the constant C > is depending only on p,9, fi, fi' , k and a. Setting 



(3.11) 



A{t) 



P 



A{t) 



P 



B{t) 



Pp{p,ey -^'■'' Pp{p,ey '-' ep,{p,e) 

thenA{t) = A{p* + a{t),'d* + 'd{t)),A{t) = A{p* +a{t),d* +d{t)) and B{t) = B{p* + a{t),d* +d{t)). 

Proof. Using the Friedrichs monifier, we may assume that (a, u;,"!?) € C(0,ii; 7^°°'°^'°°). For any 
multi- index a with < |a| < 3, applying 5" to ()3.ip i . p.ip p. p.ip ^!. then taking the Lp' inner product 
of the resultant equations with d^a{t), A{t)d^w{t), and B{t)d^i!}{t), respectively, we have 



-— \\d"a 



{ip* + ait))d^wit),Vdyit)) = {d^{v*ait)) + Ut),Vd^a{t)), 



and 



where 



(i(t)9>,(t), 5>(t)) - {Md^{^Aw{t) + {p + p')V (V • w{t))}, d^w{t)) 

+{ip* + ait))Vdy{t),d^w{t)) - K{VAdSait),A{t)dSwit)) + (i(t)i?(t)V9°^(t), 9°zx;(t)) 

= {d»fit) + Ut),Ait)d-wit)), 



{Bit)d^Mt),d^m) - a{D*Ad^m,md^m) + {Ait)B{t)Vd^m,d:iV ■ w){t)) 

= {dSh{t) + K^{t),Bit)d^m), 

Ut) = Y, Cid^^ {p* + a{t)) d^Mt), 

I3<a 

Ut) = E ^^ I f^""''^) ^-^(/^^^(*) + (^ + ^')V (V ■ w(t))) - [dT^A{t)) vd^Mt) 

/3<a ^^ ^ ^ 

-(ar^i?(t))va^^(t)}, 



I3<a 
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Canceling the terms {{p* + a{t))Vd^a{t),d^w{t)) and {A{t)B{t)Vd^i3{t),d^w{t)) by adding the 
above three formulas and using the identities 

{A{t)d^wt{t),d^wit)) = l|(i(t)9>(t),a>(i)) - ^(i,(t)9>(t),9>(i)) 

we get from integration by parts that 
d 






+ 



+ 

+ 
1 

+ 2 



{ViA{t)B{t))d^m,dSw{t)) +K{VAdSa{t),A{t)dSw{t)) 
^|(v(|^)5>(t),Va>(t))| + (^ + /.')|(v(|^).a>(t),5-(V.u;)(t)) 
+ 



(3.12) 



{dSm,A{t)dSw{t)) 
|(I„(t),V9Mi))l + 

{Mt)d^wit),d^w{t)) 



{d^h{t),B{t)d-m) 

{Ut),A{t)d-w{t))\ + \{K^it),B{t)d^m) 

{Btit)d^m,d^m) 



1 

+ 2 



Now, we estimate Ii,i = 1,2, ... ,10 term by term. First, if a = 0, employing the Hardy inequality, 
we have 



a{t) 



\Va{t)\\ <Ce||Va(t)|| 



h< \\{1 + \x\)v*\\l^ 
If 1 < |a| < 3, using integration by parts and the Sobolev inequality, we get 
/i < E ^a k^r^V • v*d^Mt) + d^^v* ■ dlVa{t),d^a{t)) 



(3.13) 



/3<« 






LS 



d^Mt) ^, \\dy 



+ 



dT^i 



L6 



9fVa(t) 



\d^.(J 



L'i 



} (3.14) 



<Ce||Va(t)||f„| 
where we have used the following inequalities (cf.[H[TU]). 

IkllLa < H'wIIls + 11^*11x6) ||w||l6 < C'l|V'u||, "iu^H^iM?) 
I2 and Is can be estimated as follows 

h <C\\{Va\Va{t),Ve*,Vm)\\L^ W'^d^mW l|9Mi)llLS 
<Ce\\Vd^m\\' 

h<c\\{i + \x\f{ya*,ve*)\\ 



\x\ 


\x\ 



+C\\iVa,Vi9){t)\\\\d>mL^\\dym 
<Ce\\{Va,Vw,Vm)\\UH,H 



L6 



(3.15) 



(3.16) 



32 



Zhengzheng Chen and Huijiang Zhao 



Shnilar to the estimate of /2, we can get 

Is <Ce||VaXt)f . (3.17) 

Now, we turn to estimate I^. We deduce from integration by parts and the equation (|3.1|) i that 



h = -K{Ad^a{t),V-[A{t)d^w{t)^ 

= -K{Ad^a{t), {VA{t)){p* + a{t))d^w{t) + A{t)V ■ {{p* + a{t))d^w{t)}) 
= -K{Ad^a{t), (Vi(t))(p* + a{t))d^w{t)) + K{Ady{t), A{t)d^at{t)) 

+K Y, C^ai^d^it), V • {d^^ip* + a{tMw{t)}) + '^(A5-a(t), V • d^iait)v*)) 

I3<a 
= h,l + h,2 + h,3 + h,4 

I^^i can be estimated as follows 

/4,i < C\\Ad-a{t)\\\\vA{t)h4{p*,<^m\L^\\d>mLe 

<C||(V/5*,Va(t),Vr,V^(t))||i||(p*,a(t))||i^||Va>(t)||||A9,"a 

< Ce ( ||VaXt)f + ||V29°a(t)"' 



(3.18) 



For 14^2, using integration by parts and the equation (j3.ip i again, we have 
h,2 - 



^±{vdy{t),A{t)vdy{t)) + ^{vdy{t),Mt)vdy{t)) 



+ (Va-a(t), VA(t)V ■ {{p* + a{t))d^w{t)}) 



+ 



^ E Cf (Va:a(t), Vi(t)V • {aX(p* + a(t))9f «;(i)}) 



I3<a 



+ {Vdyit),VA{t)V-d^{a{t)v*)) 
K d 



2dt 



{vdy{t),A{t)vdy{t)) + il, + 1| , + /|o + /t2- 



To estimate II2, we use the equation ()3.ip i and ()3.1 



7-1 

-'4,2 



-(( A,(t)ai(t) + A,(t)^i(t) ) VaXt), V9-a(t)) 



< C ||V • {(p* + ait))w{t)} + V • (^;V(t))||^.o \\Vdy 
+C \\h{t) - a{D{t) - D*)A^{t) - E{t)(V ■ w){t)\\^. 



\Vd^a 



+- 



Ka 



{Ae{t)Dit)Am^dy{t),Vd^ait)) 



<C€\\Vd;a 



+ 



Ka 



{V[Ae{t)D{t)) ■Vm.Wd'^a{t)\^) 



+Ka 



{Ae{t)D{t)V{){t) ■ Vdy{t),Ady{t)) 



<Ce\\\7d-ait)f + C\\Vm\\ 



L°° 



iva"a 



+C\\{Va*,Va{t),Vr,V^t))\\L^\\Vdy{t)\ 



Ady 

2 



<Ce\\\7d^a{t) 



(3.19) 



(3.20) 



(3.21) 
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It follows from the Holder inequality and the Sobolev inequality that 

/|,2 < C||(Va*,Va(t),V^(t),Vr)||^e \\{Va*,Vam\^e ||5>(t)||^6 ||V9> 

+C||(/5*,a(t))||^^ ||(Va*,V(T(t),V^(i),Vr)||^^ ||9°(V-ii;)(t)|| uvaser 



(3.22) 



<C76 ||v9Mi)ll +||va>(t)||- 



and 



ll2 + ll2<CY.C^\\Vd^ait)\\\\{Va*,Vait),Vm,'^^*)\\L^ 

I3<a 






L6 



+ 



L6 



8'' 



+ 



Lfi 



+ 



.a-/3 



5r''(p* + cT(t)) 



dr^a{t) ^^ d^{V-v*) 



L6 



vdSw{t) 



L^ 



(3.23) 



<C6(||vaMi)!l? + llv^(t)||f,i 

Combining ^^^-I^^M), we obtain 

Similar to the estimate of Ii, we can also get 



/4,3<Ce(^i|Va(t)||f,|+i + ||V^i)||f,| 

hA<Ce\\Va{t)\\l^^,. 
Putting (1339]) . (I3:2i]l - (13:26]) into (I3J8D gives rise to 

h < -^j^{A{t)Vdy{t),Vdy{t)) + Ce {\\Va{t)\\l^^, + ||Vz^(t)||f,| 



(3.24) 

(3.25) 
(3.26) 

(3.27) 



To estimate Iq and I7, we use Lemma 3.1. Here, we only give the detailed estimation of I7. Ip, can 
be estimated similarly. In fact, /y can be divided into the following two parts 



l7< {H^it),B{t)d^i}{t)) 

= -^7,1 + -^7,2- 

For Iy^2, using integration by parts, we have 



+ a 



[Di{t)a{t)Ad^m,md^m) 



(3.28) 



'7,2 



+ a 



{D^{t)a{t)vd^m,^ (md^m)) 



{V{D,{t)a{t)) .vd^m,B{t)d^m) 

<\\{Va\Va{t),Vr,Vm)\\L^\\ydSm\\mm\\L^ + C\\a{t)\\L^\\Vd^mf 
<Ce\\VdSmf- 
To estimate Ii^i, we use (|3.6p . If a = 0, 



(3.29) 



17,1 <c{ \\{i + \x\y(yv* ,v^a* ,ve* ,v'e* ,H,G) 



IL°= 



[1 + \x\)v*\\l. 



m 

\x\ 



ia,w,m) 



m 



|(Vo-(t), W(t), Vw{t), Vcr*, VV*, Vv*)|| , II Vw 



|v^(i)|| + lk(t)IL3||^(t)|lz.6||||v^(t)|| 

+ ||(Va*,Va(t))||||Va(t)||^3||^(t)||,.6^ 
<Ce\\{Va,Vw,S/m)\\lo,o- 



(3.30) 



34 



Zhengzheng Chen and Huijiang Zhao 



and if 1 < \a\ < 3, 



|a|+l 

i7,i<cl Yl ||(v'^+V,v'^^*,v'^+ir,v'^-i/7)||^3lk(i)ILa||5,"^(i)|| 

\a\+2 \a\+l 

|a| + l 

+ E IIv^^1Il3 II^WIIl^ \\9^m\\ + ||vi"i+'^*||^3 ll^WlUa ||a:^(t)|| 

H 



< C{e + A) ||(Va, Vw, V^)(t)||f ,_,,,, ,, + CX~^ || Vl-I^(t 



(3.31) 



Combining (|3.28p - (|3.3ip . we have 



I7<C 



e\\{Va,Vw,V^){t)\\l^^^, a = 0, 

(e + A) \\{Va, Vw, ^m)\\U+iM,H + ^'^ || V'°l^(t)||' , 1 < |a| < 3. 



(3.32) 



By using the same argument as I7, one can get 

e\\{Va,Vw,V^)it)f , a = 0, 
l6<C { „ (3.33) 

(e + A)||(Va,Vw;,V^)(i)||f,|_i,H,i„|_i + A-i||vN^(t)|| , 1 < |a| < 3. 

Similar to the estimate of Ji , it is easy to check that 



Is<Ce\\iVa,Vw,W^)it)\\ 



a\,\a\,\a\ ■ 



(3.34) 



In order to estimate /g and Iiq, we use the equations ()3.ip i and ()3.ip 9 again. In fact for Ig, 



2/q < C 



{Ap{t)at{t)d^w{t),d^w{t)) + {Ae{t)Mt)d^w{t),d^w{t)) 



I9.1 + 1< 



(3.35) 



9,1 -r -(9,2 



Using (j3.ip i ■ (|3.ip 9 and (|3.6p . Ig^i and Ig^2 can be estimated as fohows 



^9,1 = (V • {{p* + ait))w{t) + v*a{t)}, ip(t)5Xt)5>(t)) 

<C\\{Va*,Va{t),Vw{t),Vv*)\\\\{w{t),a{t),p*,v*)h,\\d^w{t)\\l, 
<Ce\\Vd°w{t)f 



(3.36) 
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'9,2 



[aD{t)A^{t) - E{t){V • w){t) + h{t), Ag{t)d^w{t)d^w{t)) 



<a 



S/D{t)Vi^{t),Aeit)d^w{t)d^w{t)) +a {D{t)V^{t),VAg{t)d^w{t)d^w{t)) 



+2a 



{D{t)V^{t),A0it)Vd^w{t)d^wit)) 



+ 
<C 



{~E{t)iV ■ w)it) + Hoit), Ae{t)d^w{t)d^w{t)) 



{\\d>m 



(3.37) 



LB 



(V'^^*, Vv*,V^a{t), Vw{t), Vi?(t), Ve*,G, H) 



X \\{a{t),w{t),d{t),v*)\\^, + ||(Vcj(t), W(t), Vcj*, W, Vi9(t), Vr, Vu;(t), Vi;*)| 
X ||(Va, V«;,V^)(t)||i6) + ||V^(t)|1^3 i|V9,"u;(t)|| ||a°u;(t)||^a } 
<C6||Va>(i)f 

Consequently, 

Ig < Ce II V5Xt)f. 

Finally, the term Iiq is estimated in way similar to that of Ig, and we have 

/io<Ce||Va^??(t)f . 



(3.38) 



(3.39) 



Combining (|3J2]) - (|3T7l) . dS^ZD, ([QSjl - dOiD . (lOHj) and ([339]) . we can obtain (|3lil and ([XTT]) . if 

we take di = min^ gNg^^^gNJ ^ p' I' '^a = '^™-(pm£S{pfi){^^{P^^)^{P^^)} ^^^ choose e and A small 
enough. This completes the proof of Lemma 3.2. 

3.3 Estimates for Va{t) and its derivatives up to V^a{t) 

Lemma 3.3. Let {a,w,'d){t) e C{0,ti;H^'^'^) be a solution to the Cauchy problem ( (g.ip -i TOj) for 
some positive constant ti. Then there exists three constants eq, Aq > and d^ > such that if e < eo 
and ||(o", tt;,'!?)(t)||4^3^3, ||(o", t(;,'!?)||jr5,5,5 < e, then it holds 



j^{w{t),Va{t)) +d3\\Va{t)f + K\\V^a{t)f <C\\{Vw,Vm)f 



and for 1 < k < 3 and any A with < A < Aq, 

^(V'=u;(t), V'=+V(t)) + ds II V'=+V(t)||' + At 



V+'ait) <C\\{Va,Vw,V^){t)\\ 



(3.40) 



(3.41) 



where the constant C > is depending only on p,6, fi, fj,',K and a. 



Proof. Using the Friedrichs mollifier, we may assume that (a, z/;,"!?) € C{0,ti;'H°°'°°'°°). For any 
multi- index a with < |a| < 3, applying d^ to ()3.ip 9. then taking the L^ inner product of the 
resultant equations with d°Va{t), we have 

{A{t)vdy{t),vdy{t)) + K\\v^dy{t)f 

= -{d^wtit), Vd^t)) + \{ds{j, ii^Awit) + (/. + fi')V (V-) w{t)]] , VdMt)) 

+ E C^\{dr^A{t)d^,Va{t)ydy{t))\ (3.42) 

+ \{d-{B{t)vm),'^dym + \{d-f{t),vdMm 

= h + h + --- + h. 
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For Ii, we deduce from integration by parts and (j3.ip i that 

h = -|(5>(t), V9:a(t)) - {d^{V-wm,d^atit)) 

= -|(a>(t), Vd^it)) + {d^iV ■ wm,d^V ■ {{p* + ait))w{t)}) 
+{d^iV-w){t),d^V-iait)v*)) 

= -^(5>(t),V9:a(t))+/i,i+/i,2 

By using the way similar to that of (j3.14p . we have 

h,i + Ii,2 < Ce \\V^dy{t)f + Ce {\\Vw{t)\\l^^, + ||Va(t)||f,|) (3.44) 

For I2, let ao < a be a multi- index with |ao| = 1, then it follows from integration by parts and the 
Cauchy inequality that 



(5r"" { ^ [f^Aw{t) + {f, + ^')V (V.) w{t)] \ , V5r"V(t)) 

|2 , ^x_li,w, ../.Ml2 



(3.45) 



h<Ce\\Va{t)\\L, (3.46) 



Using the Cauchy inequality , I3 and I4 can be estimated as follows 

/4 < A \\Vdy{t)f + CX~^ ||VT?(t)||J,| . (3.47) 

Finally, similar to the estimate of I7 in the proof of Lemma 3.2, we have 

h < C{e + A) WVdyml + CX-' IKVcT, Vw, Vi?)(t)|l J,|_i,|,|,|,|_i . (3.48) 

Combining ()3.42p - ()3.48p and summing up a with \a\ = k, we can get (|3.4ip . if we take ds = 
min^ gNg5(^ m < ''^^' > and choose e and A small enough. For a = 0, by using the same argument 
as above, we can also get (|3.40p . This completes the proof of Lemma 3.3. 

3.4 Proof of Proposition 3.2 

Let {a,w,'d){t) G C(0,ti;'H^'^'^) be a solution to the Cauchy problem (|3.ip - p.2p for some positive 
constant ti. Furthermore, we assume that ||(cr, w,i?)(i)||4,3,3, ||((T, w;,'!?)||jr5,5,5 < e, where e > is small 
enough such that we can use the results obtained in Lemmas 3.2-3.3. Set 

[a,w,^] (t) = \\a{t)f + {A{t)Va{t),Va{t)) + {1(1) wit), wit)) + (^(t)^(t),T?(t)) 

where {Ait), Ait) and Bit) are defined as in Lemma 3.2. 

Multiplying (j3.40p with a small constant Aq, then adding the resultant equation to ()3.10p . we 
have 

j^{ao[a,w,^]it) + bo{wit),Vait))} + \\iVa,Vw,Vm)\\lo,o<^ (3-49) 
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provided that e > is small enough. Here and hereafter, au,bi, > 0,u = 0, 1,--- ,3 denote some 
constants depending only on p,6, n, fi',K and a. Then summing up (|3.41|) . (|3.11|) with k = 1 and 
([339]), we get 

4 I ^ a. [W, V'w, V^^] (t) + Y, K{^"w{t), V^+Vlt)) I + |l(Va, Vw, V'd){t)\\l , , < (3.50) 
^* U=o v=Q ) 

Similarly, summing up (|3.41|) . (|3.11|) with k = 2 and (|3.5Up gives 

^ I E "- [W, V'w, V^^] (i) + Y, buiS/'wit), V^+Vlt)) I + II (Va, S/w, Vi^mWl 2 2 < (3.51) 
^^ U=o i.=0 J 

and summing up (|3.4ip . (|3.1ip with k = 3 and (|3.5ip gives 
^ I E «- [^'^' ^"^' ^'^] (*) + E &.(V^u;(t), ^^+^(0) I + II (V^, Vw, Vm)\\l3 3 < (3.52) 

"'^ U=0 i.=0 J 

Integrating (j3.52p with respect to t over [0,t], we have 

N[a,w,^]{t)+ [ \\{Va,Vw,V^){s)\\l^^ ds < N[V''a,V''w,V''^]{0) (3.53) 

Jo ' ' 

for any t G [0, ti], where 

3 3 

N [a, w, ??] (t) = E "'^ t^''^' ^''^' ^''^l (*) + E &^(V''^(t), V'^+V(t)), t > 

Denote Bq = min(^^g)g5(^ g) |i(/), 6^), i(/), 0),B{p, 6), l| and Si = max(^ ^^^^(^ g){i(/), 6^), i(p, 6^), 

-B(/9, 9), 1}. Since we may assume without loss of generality that a^, < a^-i and hy < a^ min{i?0) l}/4 
for V = 1, 2, 3, it follows from simple calculation that 



"3 „ uf .Q^/J.Ml2 / AT r .QW^^ ^ o_ \\( .a\/M\\\2 

4 



So||(a,u;,T?)(t)||^33<Ar[a,u;,^](t)<2ao||(cT,u;,79)(t)||^^3^3 (3.54) 



for each t E [0, ti]. Combining ()3.53p and ()3.54p . we get ()3.3p . This completes the proof of Proposition 
3.2. 

Hence, by Propositions 3.1 and 3.2, we finally arrive at the conclusion of Theorem 1.2. 



Acknow^ledgment 

This work was supported by a grant from the National Natural Science Foundation of China under 
contract 10925103 and "the Fundamental Research Funds for the Central Universities". 



38 Zhengzheng Chen and Huijiang Zhao 

References 

[1] R. Adams, Sobolev Spaces, Academic Press, New York, 1985. 

[2] J. E. Dunn, J. Serrin, On the thermomechanics of interstital working. Arch. Ration. Mech. Anal, 
88 (1985), 95-133. 

[3] D. M. Anderson, G. B. McFadden, G. B. Wheeler, Diffuse-interface methods in fluid mech, Ann. 
Rev. Fluid Mech., 30 (1998), 139-165. 

[4] J. W. Cahn, J. E. Milliard, Free energy of a nonuniform system, I. Interfacial free energy, J. 
Chem. Phys., 28 (1998), 258-267. 

[5] M. E. Gurtin, D. Polignone, J. Vinals, Two-phase binary fluids and immiscible fluids describled 
by an order parameter. Math. Models Methods Appl. Sci., 6 (6) (1996), 815-831. 

[6] Y. Shibata, K. Tanaka, On the steady flow of compressible viscous fluid and its stability with 
respect to initial disturbance, J. Math. Soc. Japan, 55 (3) (2003), 797-826. 

[7] J. Z. Qian, H. Yin, On the stationary solutions of the full compressible Navier-Stokes equations 
and its stability with respect to initial disturbance, J. Differential Equations, 237 (2007), 225- 
256. 

[8] D. Bresch, B. Desjardins, C. K. Lin, On some compressible fluid models: Korteweg, lubrication 
and shallow water systems, Conini. Partial Differential Equations, 28 (2003), 843-868. 

[9] B. Haspot, Existence of global weak solution for compressible fluid models of Korteweg type, J. 
Math. Fluid Mech., 13 (2011), 223-249. 

[10] B. Haspot, Existence of strong solutions for nonisothermal Korteweg system, Ann. Math. Blaise 
Pascal, 16 (2009) (2), 431-481. 

[11] R. Danchin, B. Desjardins, Existence of solutions for compressible fluid models of Korteweg 
type, Ann. Inst. Henri Poincare Anal. Nonlinear, 18 (2001), 97-133. 

[12] H. Hattori, D. Li, Solutions for two dimensional systems for materials of Korteweg type, SIAM 
J. Math. Anal, 25 (1994), 85-98. 

[13] H. Hattori, D. Li, Golobal solutions of a high dimensional systems for Korteweg materials, J. 
Math. Anal. Appl, 198 (1996), 84-97. 

[14] M. Kotschote, Strong solutions for a compressible fluid model of Korteweg type, Ann. Inst. 
Henri. Poincare Anal. Nonlinear, 25 (4) (2008), 679-696. 

[15] Y. J. Wang, Z. Tan, Optimal decay rates for the compressible fluid models of Korteweg type, J. 
Math. Anal. Appl, 379 (2011), 256-271. 

[16] Z. Tan, H. Q. Wang, X. J. Kai, Global existence and optimal L^ decay rate for the strong 
solutions to the compressible fluid models of Korteweg type, J. Math. Anal. Appl. 390 (2012), 
181-187. 

[17] Y. P. Li, Global existence and optimal decay rate for the compressible Navier-Stokes-Korteweg 
equations with external force, J. Math. Anal. Appl, 388 (2012), 1218-1232. 



Stationary solutions to the full compressible Navier-Stokes-Korteweg system 39 

[18] L. C. Evans, Partial Differential Equations, Graduate studies in Mathematics, Volume 19, 1998. 

[19] M. E. Talor, Partial Differential Equations (III), Nonlinear Equation, Springer, 1996. 

[20] A. Matsumura, T. Nishida, The initial value problem for the equations of motion of compressible 
viscous and heat-conductive fluids, Proc. Japan Acad. Ser. A Math. Sci.,55 (1979), 337-342. 

[21] A. Matsumura, T. Nishida, The initial value problems for the equations of motion of compressible 
viscous and heat-conductive gas, J. Math. Kyoto Univ., 20 (1980), 67-104. 

[22] A. Matsumura, T. Nishida, Initial-boundary value problems for the equations of motion of 
compressible viscous and heat-conductive fluids, Comm. Math. Phys., 89(4) (1983), 445-464. 

[23] A. Matsumura, T. Nishida, Exterior stationary problems for the equations of motion of com- 
pressible viscous and heat-conductive fluids, in: Differntial Equations, Xanthi,1987, in: Lecture 
Notes in Pure and Appl. Math., 118 Dekker, New York, 1989, 473-479. 

[24] A. Novotny, M. Padula, Physically reasonable solutions to steady compressible Navier-Stokes 
equations in 3D-exterior domains (t^oo = 0), J. Math. Kyoto Univ., 36 (1996), 389-422. 

[25] A. Novotny, M. Padula, Physically reasonable solutions to steady compressible Navier-Stokes 
equations in 3D-exterior domains (t^oo 7^ 0), Math. Ann., 308 (1997), 439-489. 

[26] K. Tanaka, A stability of steady flow of compressible viscous fluid with respect to intial distur- 
bance (foo / 0), Math. Methods Appl. Sci., 29 (2006) (12), 1451-1466. 



